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Foreword to the preprint version 

This report summarizes a post-doctoral training period of four months that I spent, eight years ago, in 
the group working on Propagation and Scattering of Electromagnetic Waves that is led by Jean-Jacques 
Greffet (Ecole Centrale Paris). The report deals with classical and quantum descriptions of particles 
that interact with smooth random potentials, for example ultracold atoms in the dipole potential of an 
optical speckle pattern. In addition, a discussion of the link between Radiative Transfer theory and the 
underlying coherence theory of wave optics is presented. The Radiative Transfer Equation is shown to 
emerge as a limiting case of the Bethe-Salpeter Equation, and some next-order corrections are discussed. 

I decided to put this material on the arxiv because of a recent explosion of interest in the transport 
of ultracold atoms in the multiple scattering regime. Also a number of colleagues have benefited from 
the parts dealing with Radiative Transfer. I hope that these pages may be useful for other interested 
scientists as well. 

The report is written in French which may cause difficulties, although only basic vocabulary is used, 
and the technical terms are mainly identical to the English ones (some English equivalents are given in 
footnotes, if needed). Just keep in mind that the French diffusion means 'scattering', while diffusion 
spatiale is used for 'diffusion' . In case you really get stuck, please do not hesitate to contact me at 

carsten.henkel@physik.uni-potsdam.de 

C. Henkel, Potsdam, 03 May 2005. 
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Introduction 



Le probleme 

Le present rapport resume les quelques idees que j'ai pu explorer au sujet du trans- 
port, que ce soit de lumiere ou d'atomes. Ce sujet est ne d'une question apparue vers 
la fin de ma these : 

lorsque des atomes de basse temperature traversent un champ lumineux de 
tavelures, qu'arrive-t-il a leur distribution de vitesse ? Et qu'en est-il de la 
coherence spatiale des ondes de DE Broglie associee aux atomes ? 

Une telle question a un interet immediat pour ce que Ton appelle « l'optique atom- 
ique», qu'il s'agisse de l'imagerie d'atomes ou de l'interferometrie. Le groupe d'op- 
tique atomique a l'lnstitut d'Optique a Orsay a pu observer dans l'experience du miroir 
a atomes [Landragin96b| que la traversee du champ de tavelures elargit la distribution 
des vitesses transverses, perpendiculaires a la vitesse moyenne des atomes. Dans une 
telle experience, des champs de tavelures apparaissent a cause de faisceaux lumineux 
parasites ou bien par diffusion sur une surface rugueuse. L'elargissement de la vitesse 
transverse depend en particulier de la quantite de lumiere parasite. Dans le cas le moins 
favorable, il empechait d' observer la diffraction d'atomes par un reseau en reflexion 
parce que la largeur angulaire du faisceau atomique devenait trop importante. 

D'un point de vue theorique, le mouvement d'une onde de matiere dans un champ 
de tavelures se resume a un probleme de propagation d'une onde dans un milieu alea- 
toire ou l'onde subit des multiples collisions elastiques. Sous cet elairage, Ton concoit 
aisement qu'il y a un grand nombre de situations physiques ou apparait un probleme 
analogue : 

- la diffusion de la lumiere dans un milieu fortement diffusant (atmospheres stel- 
laire ou terrestre, liquides, peintures, tissus biologiques) ; 

- le transport des electrons dans un solide « reel » (avec des defauts et des vibra- 
tions thermiques), avec des applications pour la conductivite ; 

- la diffusion des neutrons dans les reacteurs a fission ; 

- la propagation des ondes sismiques a travers le globe terrestre ; 

- et d'autres exemples encore (arret d'un faisceaux de particules radioactives dans 
une matrice solide, mouvement de boules dans un terrain rugueux, marche au 
hasard d'une particule, mouvement d'atomes dans un champ laser resonnant ...) 
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II est done naturel d' explorer les analogies avec ces problemes pour trouver des meth- 
odes qui permettent de repondre aux questions posees. 

Approches theoriques 

Une caracteristique commune de ces methodes theoriques est qu'elles sont des 
approches statistiques : au lieu de decrire le mouvement d'une particule (ou d'une 
onde) donnee dans un milieu aleatoire donnee, Ton se resigne a calculer un mouve- 
ment moyen, que Ton espere « typique » pour une mesure portant sur un grand nombre 
de particules (ondes). L'on construit ainsi une theorie de transport pour les particules 
(ondes). La quantite dont on calcule la valeur moyenne depend de la situation experi- 
mentale que Ton envisage : une fonction d'onde moyenne si celle-ci interfere avec 
un « oscillateur local », e'est-a-dire une onde de reference avec une phase fixee ; une 
fonction de coherence de l'onde si Ton etudie une figure d' interference moyennee sur 
un grand nombre de realisations du systeme ; ou encore une fonction de correlation 
d'ordre quatrieme (correlations d'intensite), lorsque Ton s'interesse par exemple aux 
proprietes de la figure d' interference creee par une seule realisation de l'experience 
(cette situation s'applique a 1' interference entre deux condensats de Bose-Einstein 
ou encore aux tavelures lumineuses d'un faisceau laser coherent). 

Les exemples que Ton vient de donner s'inscrivent essentiellement dans un cadre 
statistique ondulatoire, ou les objets de base sont des ondes (ou champs), comme les 
ondes de matiere ou lumineuses. II existe egalement des problemes ou ce sont des par- 
ticules qui sont diffusees dans un milieu aleatoire. Dans ce cas, il n'a pas de sens de 
parler d'ondes ou d' interference : la description naturelle des moyennes du systeme 
passe simplement par une densite spatiale ou une distribution des vitesses. L'on est 
habitue a une telle approche par la theorie cinetique des gaz dont le point de depart 
est l'equation de Boltzmann : cette equation decrit comment la fonction de distri- 
bution d'un ensemble de particules evolue en presence d'une force exterieure d'une 
part, et des collisions entre particules d' autre part. Un autre exemple, un peu moins 
habituel, est celui du transfert radiatif : l'on considere alors le transport de l'energie 
lumineuse dans un milieu diffusant (et/ou absorbant), en decrivant cette energie par sa 
distribution en frequence, en position et direction de propagation (la « luminance »). 
L'equation de base de cette theorie ressemble a celle de Boltzmann et s'appelle 
celle du transfert radiatif ; elle a ete introduite au debut du siecle et etudiee abon- 
damment par Chandrasekhar dans un contexte astrophysique, pour modeliser les 
atmospheres des etoiles [ Chandrasekhar |. Les ingredients de cette equation sont les 
proprietes d'absorption et de diffusion du milieu. Pour des atomes dans un potentiel 
aleatoire, une telle approche serait qualifiee de « classique », en opposition a une ap- 
proche « quantique » ou l'on manipule des ondes de matiere. Pourtant, comme on sait 
que ce sont ces dernieres qui ont une realite plus fondamentale, l'on peut se poser la 
question du rapport entre ces deux approches theoriques : comment se justifie l'ap- 
proche corpusculaire a partir de 1' approche ondulatoire? Ou encore : dans quelle lim- 



5 



ite peut-on se contenter d'une approche corpusculaire, et ou vont se manifester des 
phenomenes proprement ondulatoires ? 

Petit historique 

Ces questions ont ete poursuivies dans le contexte optique depuis les annees 
1960 environ. II s'agit de fonder l'equation du transfer! radiatif (description quasi- 
corpusculaire, absence d' interferences) sur la theorie de la diffusion multiple des on- 
des dans les milieux aleatoires. C'est maintenant que Ton se rend compte du tra- 
vail enorme des physiciens sovietiques dans ce domaine IK ravtsov96l . II est devenu 
clair que 1' approche du transfert radiatif est valable pour des champs lumineux dont 
les profils spatiaux varient lentement a l'echelle de leur longueur de coherence (des 
champs « quasi-uniformes »). Cette approche n'est plus valable, par exemple, pour 
rendre compte de la retro-diffusion exaltee ("coherent backs cattering"), ou encore de 
la localisation (« forte ») des ondes dans un milieu desordonne. Dans ce rapport, nous 
allons nous appuyer abondamment sur l'analogie avec la diffusion de la lumiere. En 
particulier, nous allons elaborer un peu les conditions de validite de l'equation du trans- 
fert radiatif. 

Dans le domaine du transport des electrons dans les solides, le developpement 
s'est poursuivi de facon plutot independante, semble-t-il. Le travail de P.-W. Ander- 
son sur la localisation des electrons dans un reseau desordonne (^'Absence of diffusion 
in certain random lattices", [ Anderson58 1) a ete d'un impact remarquable parce qu'il a 
donne un exemple ou les interferences ont des consequences dramatiques pour le trans- 
port (une propriete plutot « macroscopique »). Depuis, d'autres phenomenes d'origine 
ondulatoire ont ete etudies, comme par exemple les « fluctuations de conductance uni- 
verselles ». Un temps assez long s'est ecoule jusqu'a ce que Ton prenne conscience 
que la localisation des electrons provient simplement de leur caractere ondulatoire, 
et que la localisation existe par consequent egalement pour les ondes lumineuses (voir 
| vanTiggelen96|). C'est dans ce contexte qu'a ete predite et observee la retro-diffusion 
exaltee de la lumiere, un autre effet au-dela de l'equation standard du transfert radiatif. 

En ce qui concerne les neutrons, les recherches intenses dans les reacteurs nu- 
cleates ont rapidement conduit a une theorie du transport par analogie avec le transfert 
radiatif. Des formalismes tres divers ont ete utilises, par physiciens et mathematiciens, 
pour la justifier a partir d'une approche microscopique HOverhauser591 IOsbor n65 
Riska68 , Bobker79|. La conclusion actuelle semble que l'equation du transfert radiatif 
est valable aux grandes echelles spatiales et pour des energies pas trop basses. 

Plan du rapport 

Ce rapport contient trois chapitres dont le premier presente l'equation du transfert 
radiatif pour la lumiere d'un point de vue phenomenologique. Nous etudions en parti - 
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culier la limite ou le milieu diffuseur favorise la diffusion vers l'avant (les petits angles 
de diffusion). Dans ce regime, Ton peut reformuler le transfert radiatif a l'aide d'une 
equation de Fokker-Planck. 

Nous retrouverons cette approche au second chapitre qui est consacre au point de 
vue classique du transport d'atomes dans un potentiel aleatoire. Nous montrons que 
c'est aux echelles spatiales grandes par rapport a la longueur de correlation du poten- 
tiel, que la distribution moyenne des atomes dans l'espace des phases verifie une equa- 
tion de Fokker-Planck. Nous en calculons le tenseur de diffusion des vitesses. Aux 
temps encore plus longs (pour des atomes tres lents), la distribution spatiale atomique 
est decrite par une equation de diffusion spatiale dont nous determinons le coefficient 
de diffusion. 

Le troisieme chapitre etudie le point de vue ondulatoire : la diffusion multiple des 
ondes de matiere dans un milieu aleatoire. Nous introduisons les concepts statistiques 
pertinents ainsi que leurs equations d' evolution. Le fil rouge du chapitre est de detailler 
le passage vers une approche du type transfert radiatif. Le probleme est presente de 
facon a couvrir aussi bien la diffusion de la lumiere que celle des atomes. Nous mettons 
en evidence l'echelle spatiale minimale au-dessus de laquelle l'equation du transfert 
radiatif est equivalente aux approches ondulatoires statistiques. 
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Chapitre 1 
Transfert radiatif 



Une petite introduction dans la theorie du transfert radiatif. La lumiere y est decrite 
par la luminance. Nous etudions un peu plus en detail la situation d'une diffusion 
« piquee vers l'avant », ou l'equation du transfert radiatif peut etre ecrite sous la forme 
d'une equation de Fokker-Planck. Nous verrons au chapitre suivant qu'une equa- 
tion identique apparait pour le transport d' atonies. On illustre de cette facon le passage 
entre des theories ondulatoire et corpusculaire. 



1.1 La luminance 

La quantite centrale pour decrire la lumiere dans le transfert radiatif est la lu- 
minance Z(r, u) : elle represente la puissance rayonnee dJF par un element de sur- 
face da situe a la position r, dans un element d' angle solide df2 autour de la direc- 
tion caracterisee par le vecteur unitaire u (voir le chapitre 5.7 de Mandel et Wolf 
UMan delWolfl et le chapitre 4.3 de Rytov, Kravtsov et Tatarskii | |Rytov| ) : 



dF = J(r, u)da ■ dQ (1.1) 
L'equation d'evolution pour la luminance est celle du transfert radiatif : 

u • V r X = £(rvu) " £(r)J J d^p(u, u^; r)J(r, uQ (1.2) 

sources attenuation 4?r re-diffusion 

Cette equation traduit les processus suivants : 

- l'energie lumineuse est transportee dans la direction u le long des rayons 
geometriques ; 

- la lumiere peut etre emise par une source S(r, u) dans la direction u ; 

- elle peut etre attenuee par de 1' absorption ou de la diffusion. Ce processus est 
caracterise par un coefficient /3(r) ; 
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- la diffusion redistribue les directions de l'intensite lumineuse. Ce processus est 
decrit par le libre parcours moyen 1 £ cx (r) (egal a la section efficace totale par 
unite de volume) et la « fonction de phase » p(u, u'; r) que Ton peut interpreter 
comme une « section efficace differentielle de diffusion » pour un processus ou 
un rayon venant d'une direction u' est diffuse dans la direction u. 

1.2 Hypotheses pour le transport radiatif 

1. L'on decrit le milieu a une echelle spatiale grande devant la longueur d'onde 
optique (pour pouvoir definir une puissance rayonnee dans une direction par 
unite de surface). Ceci implique que les coefficients /3(r) et /(u, u'; r) et bien 
sur la luminance X(r, u) elle-meme varient lentement, en fonction de la position 
r, a 1' echelle de la longueur d'onde optique. 

2. La diffusion de la lumiere est supposee elastique : c'est seulement la direction de 
la lumiere, mais non pas sa frequence qui est modifiee par le milieu. Dans cette 
hypofhese, Ton peut formuler une equation du transport radiatif pour chaque 
frequence lumineuse separement. II est cependant immediat de rendre compte 
de la diffusion inelastique en introduisant une section efficace inelastique et une 
dependance de la luminance de la frequence. 

1.3 Relation a une equation de Fokker-Planck 

Afin de comparer l'equation du transport radiatif (11.21) au transport d'atomes, une 
petite modification : nous allons en deduire une equation de type Fokker-Planck 
que nous retrouverons egalement pour les atomes. A cette fin, nous supposons que la 
re-diffusion de la lumiere s'effectue surtout autour de la direction vers l'avant, c'est- 
a-dire que dans la fonction de phase /(u, u'; r), le vecteur « diffuse » u pointe dans 
une direction voisine du vecteur « incident » u'. Une telle situation correspond, par ex- 
emple, a la diffusion par des inhomogeneites d'indice, lorsque celles-ci sont beaucoup 
plus grandes que la longueur d'onde. Ce cas se presente frequemment dans la diffusion 
de la lumiere par des tissus biologiques. 

Dans la limite d'une diffusion piquee vers l'avant, nous pouvons developper la 
luminance X(r, u') au troisieme terme dans dl.21) en fonction de u'. Ceci faisant, l'on 
obtient l'expression suivante pour le terme de re-diffusion : 

re-diffusion -> — — J(r, u) + ^ F a (r)d Ua l(r, u) + ^ D a/3 (r)d Ua d u ^l(r, u) 

(1.3) 

'L'indice « ex » indique qu'il s'agit de la distance caracteristique pour l'extinction de la partie co- 
herente du champ lumineux. 
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ou u a fi designent les directions perpendiculaires au vecteur u. (Ce sont seulement 
ces directions-ci qui interviennent dans la derivee, etant donne que u est un vecteur 
unitaire.) Dans (II .3b - F a (r) est la « force moyenne » et D a p(r) le tenseur de diffusion, 
lis sont donnes par 

F a (r) = j^Jdtf(u' a -u a )p{u,u'-,T) (1.4) 

An 

F>ap{r) = (U a -U a ){u' p -Up)p{u,\i!-Y) (1.5) 



AlT 



Notons qu'une telle formulation du transfer! radiatif a ete etudiee en detail par G. C. 
Pomraning [Pomraning95 1, dans le contexte de 1' arret d'un faisceau de particules a 
l'interieur d'un solide diffusant. Le groupe de J. M. Luck l'a egalement utilisee pour 
trouver des solutions analytiques dans une geometrie planaire ILuck96l . 

Lorsque le milieu diffuseur est statistiquement isotrope (la fonction de phase ne 
depend que de cos 9 = u ■ u'), la force F a (r) s'annule et le tenseur D a p(r) est propor- 
tionnel au Kronecker 5 af3 : 

milieu isotrope : D af3 (r) = -^f-(sin 2 9) p » J^-(9 2 ) p - ^ ; 



44x(r) N " /P 44x(r) N /P 2£ t * rans (r) 

ou ^trans = 4x/ (1 — (cos 9)) est le libre parcours moyen de transport. Nous constatons 
done que les elements du tenseur de diffusion correspondent a la valeur quadratique 
moyenne de la deviation angulaire par unite de longueur des rayons. Ceci appelle a 
une interpretation physique de l'equation de Fokker-Planck (1 1 .31 ) . 



1.4 Mouvement diffusif des directions des rayons 

Introduisons les angles 9, pour le vecteur unitaire u et supposons que le milieu 
diffusant est sans sources, non-absorbant et isotrope. Le tenseur de diffusion s'ecrit 
alors Dap = D 5 a f3 et l'equation de F.-P. pour la luminance devient 



u ■ V r X(r, 9, 



D 







sin 9 89 



sin 9 



cT 
89 



+ 



8 2 



sin 2 9 , 



I(r,9, 



(1.6) 



Les derivees angulaires correspondent au Laplacien sur la sphere. Leurs fonctions pro- 
pres sont done les polynomes de Legendre P„(cos 9)e tm ^ . Pour une luminance qui 
est initialement collimatee autour de la direction 9 = (l'axe Oz), FRISCH donne 
alors la solution explicite suivante ( HFrisch 66l. eq. (19.20)) 



1 oo 
^ n n=0 



l)P„(cos( 



-n{n+l)Dz 



(1.7) 



Z(z = 0,9,^) = 5(9) 
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On en deduit que les structures angulaires fines (les harmoniques spheriques superieurs 
avec n grands) sont amorties plus rapidement que les variations angulaires « molles ». 
Dans la limite z — ► oo d'un milieu profond, il ne subsiste plus qu'une distribution 
angulaire isotrope (n — 0). La distance caracteristique L pour l'« isotropisation » de la 
luminance est donnee par 1' inverse du coefficient de diffusion 

1 /* 

L = hl^ (1.8) 

D 2 v 

Pour une fonction de phase tres piquee vers l'avant, cette distance (le libre parcours 
moyen de transport) est beaucoup plus grande que la longueur d' attenuation (~ £ ex ) 
de la partie collimatee du faisceau lumineux. La solution explicite (11.71) permet done 
d'apprecier la transition entre un « transport balistique » (distribution angulaire piquee) 
et un transport diffusif (distribution angulaire isotrope). Ce qui lui manque, e'est de 
prendre en compte le profil spatial de l'intensite incidente que Ton suppose en effet 
uniforme. 

A titre d'exemple, considerons revolution de la « directivite » de la luminance ; la 
solution (11.71) donne 

(cos0) = e - */*k»» (1.9) 
Pour des distances faibles par rapport a £* rans , on trouve done un elargissement diffusif 

(e 2 )^2z/t t ^ ns , (i.io) 

caracteristique pour une marche au hasard des directions de la lumiere sur la sphere des 
vecteurs unitaires. D'autre part, aux grandes distances z £* rans , Ton a (cos 9) = : 
la distribution angulaire est devenue isotrope. 



Remarque. L'analogie entre cette formulation du transfert radiatif et le mouvement 
d'un ensemble de particules classique a fait naitre l'idee de « simuler » la diffusion de 
la lumiere par des techniques Monte Carlo [Feld94|. L'idee est de representer la lumi- 
nance comme la somme sur un grand nombre de trajectoires qui ont subi des diffusions 
dans le milieu. Une telle approche ressemble a la formulation de la mecanique quan- 
tique en termes d'une « integrate de chemins », et Ton peut profiler de cette analogie 
pour faire un calcul analytique approche de l'integrale de trajectoires, en retenant, par 
exemple, seulement la « trajectoire la plus probable » dans la somme. Voir le papier de 
Perelman et al. [Feld94| pour davantage de references. 



1.5 Conclusions 

Dans la limite des faibles angles de diffusion (9 2 ) 1, l'equation du transfert ra- 
diatif se transforme d'une equation integro-differentielle en une equation aux derivees 
partielles. Elle ressemble a une equation de Fokker-Planck pour une distribution 
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de probabilite. Une telle description est justifiee lorsque la fonction de phase, en fonc- 
tion de Tangle de diffusion, a une portee beaucoup plus courte que la luminance en 
fonction de la direction. 

Les directions des rayons effectuent alors un mouvement diffusif sur la sphere des 
vecteurs unitaires qui est caracterise par le tenseur de diffusion D a/3 . Pour des milieux 
diffuseurs plus epais que ~ l/D, la distribution angulaire du rayonnement devient 
isotrope. 
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Chapitre 2 

Transport d'atomes — point de vue 
classique 



II s'agit maintenant de formuler le probleme du mouvement classique d'un ensem- 
ble d'atomes dans un milieu aleatoire. Avant d'entreprendre cette etude, nous allons 
preciser un modele pour le milieu aleatoire : il s'agit d'un potentiel aleatoire tel qu'il 
apparait pour le deplacement lumineux dans un champ de tavelures lumineuses. En- 
suite, nous proposons une approche classique pour le transport d'atomes : les atomes 
seront considered comme des particules ponctuelles qui se deplacent dans le potentiel 
aleatoire. Au chapitre suivant, nous utiliserons une approche ondulatoire ou Ton etudie 
la propagation des ondes de matiere dans le potentiel. Dans les deux cas, nous allons 
faire le parallele a la propagation de la lumiere. 

2.1 Modele pour un potentiel aleatoire 

Nous supposons les atomes places dans un champ de speckles lumineux. Nous al- 
lons nous restreindre a un champ lumineux monochromatique dont la frequence est 
desaccordee beaucoup par rapport a une frequence de resonance atomique. Dans cette 
situation, le champ electrique E(r) (amplitude vectorielle complexe) donne lieu a un 
potentiel lumineux proportionnel au carre du champ et inversement proportionnel au 
disaccord en frequence uo L — uoa- Nous nous limiterons dans un premier temps a un 
atome avec un seul niveau interne dans l'etat fondamental, done a une transition atom- 
ique J = — > 1. Le potentiel lumineux V(r) prend alors la forme 



il est done proportionnel a l'intensite lumineuse. Dans cette expression, d est l'element 
de matrice de l'operateur dipole electrique entre l'etat fondamental et l'etat excite. 
Une ecriture alternative du potentiel lumineux qui fait directement intervenir l'intensite 



V(r) 




h(uj L - u A ) ' 
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lumineuse I(r) = |E(r)| 2 est 

hT 2 I(r) 

V(V) = — ™ r'-P- (2.1) 

8{uj l - u A ) l sat 

ou I sat est Vintensite de saturation. On deduit de cette formule que l'ordre de grandeur 
du potentiel lumineux est (tres grossierement) donne par la largeur en energie TiT de 
l'etat excite. 

Les tavelures lumineuses ont une structure spatiale erratique dont l'echelle spatiale 
typique est la longueur d'onde lumineuse = 2n/k L = 2tcc/uj l . A cause de leur 
complexite, Ton peut tout au plus raisonner en termes de leurs proprietes moyennes ou 
statistiques. Nous allons supposer que ces proprietes sont entierement contenues dans 
\a fonction de correlation de Vintensite lumineuse : 

(V(r)V(r')) = 5 2 (/(r)/(r')> = V 2 [1 + g(\r - r'|)] (2.2) 

oil B est le coefficient de proportionnalite entre l'intensite et le potentiel lumineux, V 
est la valeur moyenne du potentiel et g(r) sa fonction de correlation (sans dimension). 
Nous supposons que les tavelures sont statistiquement homogenes (les correlations ne 
dependent que de la difference r— r' entre deux positions) et isotropes (g ne depend que 
du module du vecteur difference). L'echelle de variation caracteristique de la fonction 
de correlation est la longueur d'onde lumineuse : a des positions plus distantes que 
quelques Xl, la fonction de correlation decroit rapidement vers zero 

r > X L : g{r) -> 0. 

Un exemple particulier de fonction de correlation pour un champ de tavelures est donne 
par : 



9{r) 



sin kLr\ 2 



V k L r J 

qui correspond a un champ lumineux avec une distribution angulaire isotrope dans tout 
1' angle solide (voir GORI HGori94ll et NUSSENZVEIG | |Nussenz veig87a|). Un autre 
exemple est une fonction de correlation gaussienne : 



r 2 



9^) = e W y--j 

qui correspond au champ lumineux rayonne dans le champ lointain par une source 
incoherente non ponctuelle dont le profil spatial d'intensite est gaussien. 



Remarques. (i) Dans ce modele, le milieu aleatoire est decrit par un potentiel, a 
la difference du mouvement brownien, par exemple, ou Ton modelise les chocs que 
subit une particule immerse dans un fluide par une force aleatoire et dependante du 
temps. Dans ce probleme apparait par exemple un amortissement de la vitesse atom- 
ique, en consequence des chocs avec le fluide environnant. Dans notre probleme par 
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contre, l'energie de l'atome est conservee a priori, il peut tout au plus y avoir un trans- 
fert d'energie entre energies potentielle et cinetique, ou bien entre energie cinetique 
« orientee » (correspondant a une distribution de vitesse collimatee) et « diffuse » (cor- 
respondant a un mouvement diffus de l'ensemble d'atomes). 

(ii) Si nous nous contentons de la fonction de correlation en deux points pour 
decrire le potentiel aleatoire, ceci revient a supposer que ce dernier possede une statis- 
tique gaussienne. C'est une hypothese commode pour le calcul, mais qui peut s'averer 
insuffisante dans certaines situations experimentales (diffusion par les particules en 
suspension liquide, par exemple). V. M. Finkel'berg discute un potentiel avec une 
statistique plus generale et donne les formules necessaires pour calculer les fonctions 
de correlation [Finkelberg67|. 

2.2 Distribution dans l'espace des phases 

La quantite de base pour decrire le mouvement d'un ensemble d'atomes classiques 
est sa densite de probability dans l'espace des phases /(r, v) qui donne le nombre 
d'atomes diV se trouvant autour de la position r et ayant une vitesse v : 

dN = /(r, v) dr dv (2.3) 

Nous notons l'analogie a la luminance X(r, s) definie a l'equation (11.11) . A la difference 
de cette derniere, la distribution dans l'espace des phases depend a la fois de la direc- 
tion et du module de la vitesse atomique. Ceci traduit le fait qu'en general, les atomes 
peuvent avoir une vitesse (une energie cinetique) quelconque. En outre, lors du mou- 
vement dans le potentiel aleatoire, l'energie cinetique des atomes est modifiee. A la 
difference de la lumiere, la diffusion d'atomes n'est done pas monochromatique : la 
frequence de l'onde de DE Broglie change. Nous pouvons tout au plus esperer de 
retrouver une conservation d'energie cinetique « en moyenne » lorsque le potentiel 
aleatoire a une valeur moyenne spatialement constante. Cette loi de conservation serait 
cependant approchee et seulement valable a une echelle spatiale grande par rapport a 
l'echelle de variation des tavelures (la longueur d'onde lumineuse). 



2.3 Derivation d'une equation de Fokker-Planck 

2.3.1 Approche heuristique 

Nous partons de l'equation de Liouville pour la distribution /(r, v; t) dans l'es- 
pace des phases : 

d t f(r, v; t) + v ■ V r f + -*-F(r) • V,/ = 0, (2.4) 
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ou F(r) = — V r V(r) est la force qui agit sur l'atome. 1 Un moyen simple d'inte- 
grer 1' equation (12.41) est d'utiliser la conservation du nombre d' atonies le long des 
trajectoires classiques. Les trajectoires dans l'espace des phases (r(t), v(t)) sont les 
solutions des equations du mouvement classiques, et la conservation du nombre de 
particules s'ecrit : 

/(r(t),v(t);t) = const. (2.5) 

Cette expression permet de calculer la fonction de distribution a des instants ulterieurs 
si Ton se la donne a un instant initial. 

Utilisons maintenant 1' equation (12.51) pour deux instants de temps voisins t et t — 
At : 

/(r , v ; t) = /(r(t - At), v(t - At); t - At) (2.6) 

de sorte qu'entre t — At et t, les trajectoires classiques sont de la forme (avec < t < 
At) 

r ~ r(t - t) + rv 

v ~ v(t-T) + ^fdr'F[r(f-T')] (2.7) 
M Jo 

Nous supposons done l'intervalle At suffisamment court et le potentiel aleatoire suff- 
isamment faible pour que les trajectoires ne soient pas fortement courbees sous l'in- 
fluence du potentiel aleatoire. En admettant en outre que la fonction de distribution 
/(r, v; t) varie lentement a l'echelle des deplacements (12.71) . nous pouvons la develop- 
per a 1' equation (12.61) . pour trouver 

/(r , v ; t) - /(r , v ; t - At) = 

-Atv - V r /(t - At) - (v(t - At) - v ) ■ V v f(t - At) + 

+ ~ E( v (^ " At ) ~ v o),(v(t - At) - v )A<V(t - At). (2.8) 

Nous allons maintenant prendre la valeur moyenne de cette equation par rapport au 
potentiel aleatoire (qui apparait dans la vitesse v(t — At)). Pour factoriser la valeur 
moyenne de la troisieme ligne de (I2.8I) . nous supposons que le long des trajectoires, 
la force aleatoire et la fonction de distribution forment un processus markovien : la 
fonction de distribution a l'instant t — At n'est alors pas correlee avec les valeurs de la 
force pour t — r > t— At. Cette hypothese semble raisonnable si pendant l'intervalle de 
temps At, l'atome traverse un grand nombre de longueurs de correlation du potentiel : 

vAt > \ L (2.9) 

II s'ensuit que cette approche ne peut ni decrire des atomes tres lents, ni rendre 
compte de leur mouvement a une echelle de temps comparable au temps entre deux 
« collisions ». 

'En ecrivant J2.4L nous negligeons les interactions entre atomes. 
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Avec l'hypothese de Markov, la moyenne sur les realisations du milieu aleatoire 
fait apparaitre 1' integrate de la fonction de correlation de la force dans le dernier terme 
de l'equation (12.81) . Celui-ci devient alors 

At 

— J dr[dr' 2 £ (F,[r - r^VojF^ro - r^v ]) d Vi d Vj (f) (2.10) 
o ^ 

Par consequent, le coefficient de diffusion de la vitesse est donne par l'expression 

1 A * 

A 7 (v)= lim — [dridrLGiMri - r')v] (2.11) 

ou T ^> \l/v est l'echelle de temps a laquelle Ton regarde la fonction de distribution 
moyenne (beaucoup plus grande que le temps entre deux collisions), et ou Gij(r) est 
le tenseur de correlations de la force aleatoire 

G y (r) = ^(0)^)) . (2.12) 

Avant de calculer en detail le coefficient de diffusion de la vitesse Dy(v), don- 
nons l'equation de Fokker-Planck que Ton trouve par cette procedure. II y appa- 
rait une complication parce que le coefficient de diffusion depend de la vitesse. Nous 
suivons 1' article de Hodapp et al. IDalibard95l (qui suit le livre de VAN Kampen 
| vanKampenl) pour la forme de l'equation de F.-P. : 

d t (f(r, v;f)> + v • V r (/> = {"^</> + Ai(v)^(/>| + ^^(/) 

(2.13) 

ou la force de derive Ff(v) ainsi que le dernier terme proviennent du re-arrangement 
des termes sous les derivees d Vi . La force de derive est donnee par 

Ff{v) = M y £9v j D ij (-v). (2.14) 



Calcul du coefficient de diffusion. 

de la force aleatoire (12.121) : 



Nous donnons d'abord la fonction de correlation 



ldg 

r dr 



d 2 g TiTj 



r = r 



(2.15) 



resultat que Ton deduit de la definition elementaire de la force comme le gradient du 
potentiel, et en echangeant l'ordre de la moyenne statistique et du processus limite pour 
le gradient. Ce calcul suppose que le potentiel aleatoire admet une derivee, hypothese 
qui semble justifiee pour un champ de speckles. 

Pour calculer le coefficient de diffusion ( 12.1 II ). passons aux variables d'integration 
[t[ +T2) /2, t' 2 — t[ . Pour l'integration sur la difference, nous utilisons l'equation (12.91 ) : 
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l'intervalle At est beaucoup plus long que le temps de correlation du tenseur de cor- 
relations Gij[(r[ — t^v]. Par consequent, celui-ci s'annule aux bornes d'integration 
±(t{ + T2V2 que nous pouvons done remplacer par ±00. L'integrande ne depend 
alors plus de la somme, et l'integrale est proportionnelle a At. Nous trouvons ainsi le 
resultat : 

+00 

D »W = 2^/ drG ««™^ = v(*«-^) <2 ' 16) 



-00 

00 



V 2 7dr ( dq\ 

avec K = — / — — - > 0. (2.17) 



M 2 J r \ dr 

v ' 

Le deuxieme terme dans (12.151) (qui fait intervenir la deuxieme derivee de la fonction 
de correlation) ne contribue pas au coefficient de diffusion parce que son integrate radi- 
ale s'annule. Pour les deux fonctions de correlations proposees ci-dessus, la constante 
K prend les valeurs suivantes : 

lumiere isotrope : K = —^7^7, (2.18) 

3 M l 

fWl V 2 

speckles gaussiens : K = ^YJp' <y2A9 ^ > 

Force de friction. Du tenseur de diffusion (12.161) nous deduisons que la force de 
derive F d (v) dans 1' equation de F.-P. (12.131) est en fait une « force de friction » avec 

F d (v) = -2MK^r, V v -F d = 8irMK5(v) (2.20) 

v 6 

Puisque notre approche n'est pas valable pour des atonies a vitesse nulle [voir la con- 
dition (12.91) 1. nous allons negliger le terme proportionnel a V„ ■ F d oc 5(~v)(f) dans 
l'equation de F.-P. d2~T3l . 

2.3.2 Discussion physique 

Amortissement de la vitesse moyenne. On deduit facilement de l'equation de 
Fokker-Planck (12.131) l'equation d'evolution pour la vitesse moyenne (apres deux 
integrations par parties et en negligeant les termes de bord) : 



d(v k ) 
dt 



d f 

— J drdv v k (f(r,v;t)) 

(/(r,v;t)) (2.21) 



drdv 



/ ' : ' : '' (V) ' T,d Vj D jh (v) 



M 

j 



-4K |drdv^(/(r,v;t)) (2.22) 
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La valeur moyenne de la vitesse diminue done sous l'effet de la force de friction et 
de la diffusion. Ceci se produit a une echelle de temps caracteristique T ~ v 3 /K, 
qui depend du module de la vitesse initiale. Pour que notre theorie soit valable, il faut 
que ce temps soit beaucoup plus grand que l'intervalle entre deux collisions. Cette 
condition se traduit par l'inegalite : 



II faut done que le potentiel aleatoire soit en moyenne beaucoup plus faible que l'en- 
ergie cinetique des atomes. Dans ce regime, l'amortissement de la vitesse est lente a 
1' echelle des collisions avec les speckles. 

Elargissement diffusif des composantes de vitesse transverses. Nous nous at- 
tendons a ce que l'amortissement de la vitesse moyenne s'accompagne d'une 
« isotropisation » de la distribution angulaire de la vitesse parce qu'il semble peu prob- 
able que la distribution puisse etre « refroidie » par le potentiel aleatoire. 

Nous rappelons que le coefficient de diffusion dans l'equation de F.-P. donne la 
valeur quadratique moyenne du changement de la vitesse : 



En utilisant le tenseur de diffusion (I2.16I) . nous constatons que les composantes paral- 
lels a la vitesse initiale v(0) ne subissent aucune diffusion 2 (Avn(t)} = 0. Par contre, 
les composantes perpendiculaires a v(0) s'elargissent de facon diffusive 



La distribution des vitesses devient isotrope lorsque les composantes transverses sont 
du meme ordre de grandeur que la vitesse initiale. Nous trouvons que ce processus 
a lieu sur la meme echelle de temps T ~ K/v 3 (0) que l'amortissement de la valeur 
moyenne de la vitesse. 

Variation de l'energie cinetique moyenne. Calculons maintenant une equation 
d'evolution pour la valeur moyenne de l'energie cinetique des atomes. L'equation de 
F.-P. (12.131) donne, apres un calcul similaire a celui de req. (l2.21b . 




(2.23) 



{Avi(t)A Vj (t)) ~2Aj[v(0)]t. 



(2.24) 




(2.25) 




= |drdv vF d (v) + M^^(Ai(vK) (/(r,v;*)> (2.26) 




2 Precisons que Auy denote la partie fluctuante de la vitesse. C'est la vitesse moyenne qui est amortie 
par la force de friction F d (v). 
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L'energie cinetique moyenne diminue done sous l'influence de la force de friction 
[le premier terme en crochets de (12.261) 1. Ceci apparait egalement dans l'equation du 
mouvement pour la vitesse sous 1'influence de la friction : 



dt 



F d (v) 



-2K- 



v(f) 



v(o) r 

v(0) 



v 3 (0) - 6Kt 



1/3 



(2.28) 

(0 < t < v 3 (0)/QK^j (2.29) 



Toutes les vitesses sont done amorties radialement et deviennent nulles en un temps fini 
T v (p) = v 3 (0)/6K de l'ordre de l'echelle de temps caracteristique T trouvee ci-dessus. 

Cependant, nous trouvons ici des resultats en contradiction avec 1' intuition que le 
potentiel aleatoire conserve en moyenne l'energie cinetique. Cette contradiction nous 
amene a la remarque suivante : 



Demonstration « malhonnete ». D'apres Keller (cite dans |Frisch661), une 
telle demonstration de l'equation de F.-P. pour la fonction de distribution 
moyenne est « malhonnete ». Nous faisons en effet des hypotheses sur la nature 
markovienne de la force aleatoire qui ne sont justifiees que dans un certain regime de 
parametres ou a certaines echelles de vitesse et de position (grande devant la longueur 
de correlation). Nous verrons au prochain paragraphe que Ton peut trouver l'equation 
de F.-P. par une procedure differente et mieux justifiee d'apres Frisch. Elle aura 
en outre la vertu de donner une equation de F.-P. « physiquement acceptable » qui 
conserve l'energie cinetique moyenne. 



2.4 Demonstrations plus rigoureuses 



2.4.1 Methode pedestre 

La facon la plus simple d'eviter la contradiction du refroidissement dans un po- 
tentiel aleatoire consiste a renverser une partie de 1' argumentation precedente pour de- 
duire l'equation de Fokker-Planck. (C'est ainsi que Frisch la trouve habituelle- 
ment.) Definissons d'abord les valeurs moyennes 



Ff(v) 



M lim 

At^O 



(Vi(t + At)-Vi(t)) 



At 



lim 

At-*0 



v(t) 

(Avi(t + At)Avj(t + At)) 



L'equation de F.-P. est alors donnee par 

d t (f(r,v;t))+v-V r (f) = J2d Vi 



2At 



M 



v(t) 



(2.30) 
(2.31) 



(2.32) 
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ou la force de derive F d (v) et le tenseur de diffusion -Dij(v) sont donnes par (12.301) 
et dCTl) . 

Pour calculer la derive et le tenseur de diffusion dans un potentiel aleatoire faible, 
nous pouvons integrer le gradient du potentiel le long d'une trajectoire droite, comme 
a l'equation (12.71) . De cette facon, Ton trouve que 

- la force de derive Ff(v) s'annule (en moyenne, le potentiel aleatoire ne donne 
pas de transfert de vitesse), et 

- un resultat identique a (12.161) pour le tenseur de diffusion. 

II reste cependant vrai que la valeur moyenne de la vitesse est amortie : l'equa- 
tion (12.211) est encore correcte et la seule difference par rapport au resultat (12.221) est 
un facteur numerique qui change (remplacer AK par 2K). Et finalement, pour l'en- 
ergie cinetique moyenne, l'equation (12.261) reste vraie (elle ne depend que de la forme 
de l'equation de F.-P.) et donne un resultat nul parce que le tenseur de diffusion est 
« orthogonal » a la vitesse, £ • Dij(v)vj = 0. Nous obtenons done le resultat attendu 
que I'energie cinetique moyenne n'est pas modifiee par le potentiel aleatoire. Ce qui 
se produit e'est bien une redistribution de la vitesse d'une composante collimatee vers 
une composante diffuse. 

2.4.2 Developpement multi-echelles de Ryzhik, Papanicolaou et 
Keller 

Idee du developpement. Nous allons presenter ici une methode alternative pour 
etablir l'equation de Fokker-Planck. Elle met en relief le fait que le potentiel alea- 
toire V(r) d'une part, et la fonction de distribution moyenne (/(r, v; t)) d'autre part, 
varient sur des echelles spatiales tres differentes. II s'agit alors d'eliminer les quantites 
avec des variations spatiales « rapides » et de formuler une equation fermee pour les 
quantites moyennes avec une variation spatiale « lente ». Ryzhik, Papanicolaou et 
Keller [ Keller96 1 introduisent a cet effet le rapport entre les echelles caracteristiques 
« rapides » et « lentes », ce qui revient pour notre situation a introduire le parametre 




(2.33) 



lis considerent ensuite un potentiel aleatoire faible proportionnel a y/e. Le calcul con- 
siste a effectuer un developpement asymptotique pour e«l. 

Cette methode semble d'abord se limiter a une situation ou existe une relation entre 
la force du potentiel, d'une part, et la portee des correlations, d'autre part. Cependant, 
Ton peut aussi la voir sous un autre point de vue et choisir independamment la force du 
potentiel et la longueur de correlation \ L : I'approche de Ryzhik fixe alors Vechelle 
L a laquelle on regarde la fonction de distribution moyenne. La discussion precedente 
confirme que ce point de vue est adapte a notre probleme : plus precisement, nous 
pouvons estimer que l'echelle spatiale caracteristique L de la fonction de distribution 
moyenne est egale a la distance parcourue necessaire pour que la diffusion ait rendu 
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isotrope la distribution angulaire de la vitesse 



v 4 (Mv 2 \ 
L^T~-~A L (_J (2.34) 

Cet argument permet d'etablir le lien suivant entre les petits parametres de notre prob- 
leme : 

( v V 

(2.35) 



\Mv 2 J ' 

et la condition e <C 1 est bien celle d'un potentiel aleatoire faible que nous avons ren- 
contree et utilisee ci-dessus. Lorsque Ryzhik et al. parlent d'un « potentiel aleatoire 
proportionnel a y 7 ? », ce nombre donne done l'ordre de grandeur de V(r) en unites de 
l'energie cinetique des atonies (la seule echelle d'energie dans le probleme). 

Notons que le developpement multi-echelles est a manipuler avec une certaine pre- 
caution : il faut identifier d'avance les echelles caracteristiques spatiales ainsi que les 
energies typiques du probleme et regrouper ces parametres dans une seule petite quan- 
tite e. 

Developpement. En fonction du petit parametre e, la force aleatoire est proportion- 
nelle a 

F ~ — oc y~ = —= (2.36) 



Par consequent, Ryzhik et al. ecriraient l'equation de Liouville (12.41) sous la forme 
d t f(v, v; t) + v ■ V r / + — ^F • V,/ = 0, (2.37) 

ou F est une quantite « de l'ordre unite ». La fonction de distribution est alors develop- 
pee en une serie de puissances 

/(r,v;f) = / (0) (r,v;t) + 

+ ^f {1) (r, p, v; t) + efW ( r , p,v;t) + . . . (2.38) 

ou r, p denotent les variables spatiales « lentes » et «rapides ». On fait l'hypothese 
que le premier terme /(°) est independant de la variable « rapide » p et que e'est le 
seul a donner un resultat non nul apres moyenne statistique sur le potentiel. Les termes 
suivants p l,2 > dependent a la fois de la variable rapide p et de la variable « lente » r. 

La difference en echelle entre ces deux variables se traduit par l'ecriture suivante 
de l'operateur gradient spatial dans (12.371) 

V r i-> V r + -V p (2.39) 

e 
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Nous pouvons maintenant inserer le developpement de la fonction de distribu- 
tion (12.381) dans 1' equation de LlOUVlLLE (12.371) et separer les differents ordres en 
e. Les deux ordres les plus bas donnent alors : 

ordre 0(1/ y/e) : = v ■ V p / (1) + F ■ V,J (0) (2.40) 
ordreO(l): = (d t + v • V r ) / (0) + v ■ V p / (2) + F ■ V v f {1) (2.41) 

La premiere equation permet d'exprimer la « partie fluctuante » /W de la distribution 
en fonction de la force aleatoire et de la distribution « lente » f(°\ La deuxieme equa- 
tion permet de trouver une equation fermee pour la distribution moyenne (/) = (f^) 
en inserant la solution de la premiere equation dans la valeur moyenne de la deuxieme 
equation : le terme en disparait en effet apres la moyenne. 

L'equation (12.401) peut se resoudre par une transformation de Fourier, et Ton 
trouve 

P\r,K,v;t) = — /dpexp(-iK-p)F(p) • V„/^(r,v;t), (2.42) 

IV ■ K, + T] J 

ou i] — >• + est un parametre de regularisation. Nous inserons ce resultat en (12.411) et 
trouvons apres la moyenne une equation de Fokker-Planck 

(dt + v V r ) (f) = ^^Ai(v)^(/), (2.43) 



it 



ou le tenseur de diffusion est donne par 



Comme le potentiel aleatoire est statistiquement homogene, le tenseur de diffu- 
sion (12.441) ne depend pas de p. En outre, l'integrale sur n peut se calculer 3 et Ton 
trouve [en utilisant (12.121) 1 

1 ? 

— oo 

L'on trouve alors un resultat identique a (12.161) si la fonction de correlation Gjj(r) ne 
depend pas du signe de r (c'est l'invariance par parite). En plus, le tenseur de diffusion 
se trouve « a la bonne position » dans l'equation de F.-P. ( 12.431) (entre les derivees d Vj ) 
de sorte que celle-ci est compatible avec la conservation en moyenne de l'energie. 

3 Le resultat est 

d« exp(-iK • r) _ J 6(Rj_)/v si r • v < 0, 



(27r) 3 iv ■ k + T) [0 sinon. 
ou Rj^ denote les composantes de r perpendiculaires au vecteur v. 
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Remarque. Nous constatons que le developpement multi-echelles conduit de facon 
generique a une equation de F.-P. a partir de l'equation de LlOUVlLLE. Ce type de 
description du transport d' atonies dans un potentiel aleatoire apparait done comme 
l'outil standard dans le point de vue classique. 



2.4.3 Raisonnement formel de Frisch 

Idee. Cette approche procede par analogie avec les equations de propagation des 
ondes en milieu aleatoire, en les formulant d'une facon assez generale. L'on identifie 
des « operateurs de LlOUVlLLE » £ et L\ dont le deuxieme contient la partie aleatoire 
de l'equation de propagation : 

d t A(t) = C A(t) + d(t)A(t), (2.46) 

ou A(t) designe un « champ » quelconque. Dans l'equation de LlOUVlLLE proprement 
dite (12.4b . nous avons les correspondances suivantes 

A(t) <— > /(r,v;t) 
C < — > —v ■ V r 

d(t) <— > ~F.V„ (2.47) 
Pour l'equation de propagation des ondes de SchrOdinger, l'on aurait par exemple 

A(t) <— ► ijj{r-t) 

d(t) > ~V(r) (2.48) 

Frisch montre alors que l'on peut ecrire des equations separees pour la valeur 
moyenne (A(t)) et la partie fluctuante 5A(t) = A{t) — (A(t)) du champ : 

d t (A(t)) = C (A(t)) + (d(t)6A(t)) n , q . 
d t 5A(t) = £ MW + 5{A(t)[(A(t)) + M(t)]}, 

La deuxieme equation peut etre resolue par un developpement iteratif, au moins de 
facon formelle. Sa solution exprime done la partie fluctuante du champ en fonction du 
champ moyen et du Liouvillien aleatoire. En la reportant dans la premiere equation, 
l'on trouve une equation pour le champ moyen. Dans le cas de la propagation des 
ondes, l'equation resultante est l'equation de Dyson pour le champ moyen (ip). 



Approximation pour le champ fluctuant. II est cependant difficile de re- 
soudre explicitement pour la partie fluctuante du champ. Ce que Frisch appelle 
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1' approximation de regularisation au premier ordre revient alors a ne retenir que le 
premier terme de 1' iteration. 4 L'on a alors la solution 

t 

SA(t) ~ J dt' exp [£ (t - 0] {A(0W>1 • ( 2 - 5 °) 

to 

Dans l'approche de Ryzhik et ah, une telle approximation arrive naturellement parce 
que la partie fluctuante est calculee par un developpement perturbatif : a 1' ordre le plus 
bas, c'est le Liouvillien fluctuant et le champ moyen (d'ordre zero) qui determinent les 
fluctuations du champ (voir 1' equation (12.401) ). 

Avec cette approximation, l'equation d'evolution pour le champ moyen devient 

t 

d t (A(t)) = C (A(t)) + J dt' (£ a (t) exp [£ (t - 0] A(0) W)>- (2-5D 

to 

ou nous voyons apparaitre sous 1' integrate la fonction de correlation du Liouvil- 
lien fluctuant, prise pour les instants t et t'. Entre ces instants, l'exponentielle 
exp [£ (t — £')] decrit le mouvement « libre» du systeme sous la seule influence du 
Liouvillien regulier. 

Application a notre probleme. En utilisant les correspondances (12.471) . le deuxieme 
terme a l'equation ( 12.511 ) pour le champ moyen s'ecrit 

t 

d Vi J dt' (FitfFjir + [t' - t)v)> d v . (/(r + (f - t)v, v; if)). (2.52) 

to 

ou r + {t' — t)w decrit le « vol libre» aux instants t' anterieurs a t. Nous consta- 
tons que l'equation pour la fonction de distribution moyenne n'est pas equivalente 
a l'equation de F.-P. (12.431) parce qu'elle fait intervenir une integrale sur (/) a des 
instants anterieurs. C'est seulement lorsque nous invoquons une hypothese supple- 
mental que nous retrouvons la theorie precedente : cette 1' hypothese etant que la 
portee de la fonction de correlation de la force aleatoire (en fonction de (£' — £)v) est 
beaucoup plus courte que l'echelle spatiale caracteristique pour la fonction de distribu- 
tion moyenne. Dans cette situation (que nous pourrions appeler « markovienne ») l'on 
peut sortir la fonction de distribution moyenne de l'integrale sur t' . En rempla§ant en- 
suite la borne inferieure t de l'integrale par — oo, l'on trouve l'equation de Fokker- 
Planck g!35b . 

Notons que l'approche de Frisch donne elle aussi un tenseur de diffusion a la 
bonne position (entre les derivees d Vi ). 

4 Notons que dans le contexte de la propagation des ondes, cette approximation est celle de Bour- 
RET ou encore la ladder approximation pour l'operateur de masse. Dans d'autres contextes, on l'ap- 
pelerait « theorie effective a basse energie ou a grande echelle » (physique des particules) ou encore 
« elimination adiabatique des variables rapides » (atome pilote par un champ a grand desaccord : reduc- 
tion des equations de Bloch optiques a la population de l'etat fondamental...) 
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2.5 Conclusions 



Dans le point de vue classique, le transport d'atomes dans un potentiel aleatoire 
faible est decrit de facon naturelle par une equation de Fokker-Planck, et non pas 
par une equation de transfer! radiatif. L' equation de F.-P. fait seulement intervenir la 
fonction de distribution au meme instant (elle est locale dans le temps) lorsque les 
fluctuations du potentiel aleatoire ont une echelle spatiale beaucoup plus courte que la 
fonction de distribution moyenne. Le transport est alors decrit par une equation lineaire 
aux derivees partielles. Son tenseur de diffusion est relie au tenseur de correlation de 
la force aleatoire, il depend de la vitesse des atonies. L' equation de F.-P permet de 
montrer que dans le potentiel aleatoire, la valeur moyenne de la vitesse decroit parce 
que la distribution angulaire des vitesses devient isotrope. Le temps caracteristique 
pour ce processus diffusif est plus long que le temps necessaire pour traverser une 
longueur de correlation du potentiel aleatoire. A cette echelle de temps, la distribution 
angulaire d'un faisceau d'atomes initialement bien collimate, s'elargit pour devenir 
isotrope, tout en conservant son energie cinetique moyenne. 

Nous avons constate que le tenseur de diffusion depend de la vitesse. II faut alors 
faire attention a sa position dans l'equation de F.-P. II y a plusieurs methodes pour 
trouver la « bonne equation » : un retour a la signification statistique de l'equation 
de Fokker-Planck est deja fructueux (notre « methode pedestre ») ; des approches 
plus formelles sont soit le developpement a plusieurs echelles (methode de Ryzhik, 
Papanicolaou et Keller), soit l'analogie a 1' approximation de Bourret de la 
theorie de propagation des ondes en milieu aleatoire (methode de Frisch). 

Finalement, une remarque sur les approches du type F.-P. et transfer! radiatif : 
nous pouvons attribuer leur difference aux descriptions classique ou ondulatoire du 
mouvement. En effet, la theorie classique part de l'equation de LlOUVlLLE ou intervi- 
ennent des derivees de la fonction de distribution dans l'espace des phases. L'equation 
de F.-P. herite alors de cette meme structure. Dans l'approche ondulatoire, la diffu- 
sion par le potentiel aleatoire est decrite par des integrates de diffusion (voir l'equation 
du transfert radiatif (11.21) ') parce que pour la diffusion d'une onde, tout potentiel, aussi 
faible soit-il, peut donner lieu a des grands angles de diffusion avec une probabilite non 
nulle (bien que faible). Dans une approche classique, par contre, la vitesse change tou- 
jours de facon continue, et ce changement tend vers zero lorsque le potentiel est faible. 
Les deux descriptions se rejoignent si la longueur d'onde est beaucoup plus petite que 
la taille des inhomogeneites du potentiel : la diffusion est alors piquee vers l'avant, et 
Ton peut passer du transfert radiatif a l'equation de F.-P. (voir au chapitre[T]). 

2.6 Exemples de transport classique 

Pour finir ce chapitre sur le transport classique d'atomes, nous donnons quel- 
ques exemples ou Ton peut trouver des solutions explicites a l'equation de Fokker- 
PLANCK. II s'agit d'une part de l'elargissement angulaire de la distribution des vitesses 
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Fig. 2.1 - Transport d'un faisceau d'atomes a travers une tranche de tavelures. Illus- 
tration de la geometric 

et d' autre part de la diffusion spatiale. 

Notre point de depart est l'equation de F.-P suivante 

(^ + v-V r )/(r,v;t) = ^A( T )/ (2.53) 

ou la fonction de distribution moyennee dans l'espace des phases est notee /, pour 
simplifier, et ou A[ T ) est l'operateur differentiel Laplacien transverse. En fonction des 
angles 9, 4> du vecteur v, il est donne par l'expression en crochets de l'equation (11.61) . 
La constante K dans (12.531) est proportionnelle a V 2 , elle est donnee par (12.171) . 



2.6.1 Diffusion des vitesses 

Etudions d'abord revolution d'un jet d'atomes collimate qui traverse un champ de 
tavelures. Dans un premier temps, supposons que les tavalures occupent une tranche 
d'epaisseur Z et que le jet a un profil spatial uniforme (voir la figure l2~TT) . En nous 
placant dans le regime stationnaire, la fonction de distribution ne depend que de la 
distance z et de la vitesse v. Nous notons egalement que dans l'equation de Fokker- 
Planck (12.531) . le module de la vitesse v intervient seulement comme parametre. 
Nous pouvons done raisonner pour une classe de vitesses donnee. A l'entree du champ 
de speckles, la fonction de distribution est de la forme 

f (v) 00 

f(z = 0; v) = f i (v)S(9) = i±lj2(2n + l)P n (cos9) (2.54) 

47F n=0 

ou P n (cos 9) sont les polynomes de Legendre qui permettent de representer la fonc- 
tion 5 angulaire (voir aussi (11.71) ; 9 denote Tangle par rapport a la direction initiale du 
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faisceau). En introduisant les coefficients f n (z) pour le developpement en harmoniques 
spheriques de la fonction de distribution, nous trouvons l'equation suivante : 



£cos0^(2n + l)P n (cos0) = -^£n(n + l)/ n (z)(2n + l)P n (cos0) (2.55) 



Lorsque Ton projette cette equation sur un polynome de Legendre P m (cos6>), 
le cos 9 au membre gauche introduit un couplage entre les polynomes voisins 
P m ± 1 (cos9). L'on trouve alors un systeme lineaire infini d'equations differentielles 
ordinaires que Ton peut resoudre, par exemple, avec des techniques de diagonalisation 
de matrices. Nous allons nous contenter d'une discussion qualitative du regime ou le 
champ de speckles a une epaisseur faible. Plus precisement, nous supposons que la 
distribution angulaire reste assez bien collimatee de sorte que nous pouvons remplacer 
le cos 9 dans (12.551) par 1' unite. Dans cette approximation, les differentes harmoniques 
spheriques se decouplent, et Ton retrouve la solution (11.71) donnee par Frisch 

f n (Z)=e- Kn( - n+1 ^ v4 f i (v) (2.56) 

Comme nous l'avons vu pour le transfer! radiatif pique vers l'avant (chap.d), les har- 
moniques eleves (les structures angulaires A9 fines) disparaissent sur une epaisseur 
caracteristique donnee par 

I KZ A „ (KZ\ l/2 

~ 1 Afl ~ I — ;- ) (2.57) 

A9 2 v 4 V v A J 

La solution (12.561) est valable tant que cette largeur angulaire reste faible. L'equa- 
tion (12.571) permet de traduire ce critere par une epaisseur maximale du milieu alea- 
toire : celle-ci est evidemment beaucoup plus grande que la longueur de correlation (la 
longueur d'onde optique). 

Nous passons maintenant a un deuxieme exemple ou le potentiel lumineux est 
« branche » pendant un temps T fini. Prenons encore un jet d'atomes collimate a l'in- 
stant t — et spatialement uniforme. Nous supposons egalement que le potentiel 
lumineux (lorsqu'il est branche) remplit tout l'espace. La fonction de distribution ne 
depend alors que de t et v. Nous introduisons les coefficients f n (t) d'un developpe- 
ment en polynomes de Legendre, et (12.551) devient 

J2^(2n+ l)P n (cos 0) = -4Z < n + l )Uz){2n + l)P n (cos 9) (2.58) 
n dt v n 

Dans cette equation, les harmoniques spheriques se decouplent exactement, et nous 
trouvons une solution formellement identique a (12.561) : 

f n (T) = e- Kn ^ T / v3 f l (v) (2.59) 

De nouveau, nous identifions un temps d'interaction caracteristique T iso ~ v 3 /K au 
bout duquel la distribution angulaire devient isotrope. Si Ton sait resoudre des details 



28 



angulaires A6 plus fins, l'effet du potentiel aleatoire se fait voir deja a des temps plus 
courts, de l'ordre de A6 2 x T iso . 

Notons finalement que le temps caracteristique T iso depend a la fois de la vitesse 
des atomes et du potentiel aleatoire (qui apparait dans la constante K). L'on peut done 
faire varier ces deux parametres pour explorer une large gamme de temps d'interaction 
« effectifs », meme si pour des raisons experimentales (effet de la gravite, taille de la 
chambre a vide), le temps reel est limite. 



2.6.2 Diffusion spatiale 

En derniere application du transport d'atomes, nous allons etudier le regime des 
temps tres longs, ou la distribution des vitesses est deja devenue isotrope. Les atomes 
continuent cependant d'evoluer a cause d'un gradient spatial dans la densite atomique. 
Dans ce regime, la distribution spatiale F(r; t) des atomes verifie une equation de la 
diffusion que nous allons determiner maintenant. 

A cet effet, nous ecrivons la distribution dans l'espace des phases comme la somme 
d'une partie isotrope F(r; t) et d'un terme « dipolaire » v • j(r; t) 

/(r,v;t) = F(r;t)+v-j(r;t) (2.60) 

(Cette approximation est appelee « Pi »parce qu'elle se contente de Po,i(cos#) dans 
le developpement en harmoniques spheriques.) En projetant l'equation de Fokker- 
PLANCK sur les harmoniques spheriques d'ordre n = 0, 1, Ton trouve 

d t F(r;t) + J^V r -v[vi(r;t)} = (2.61) 

2K 

d t j{r;t)+V r F(r;t) = — jr;f (2.62) 

v 

En nous pla§ant a des temps longs par rapport au temps caracteristique d' evolution T iso 
de la distribution angulaire, nous pouvons nous contenter de la solution stationnaire de 
l'equation (12.621) : elle relie le « flux diffusif » j(r; t) au gradient de la densite spatiale 
(loi de FlCK) 

j(r;t) = -^V r F(r;t) (2.63) 

En reportant ce resultat a l'equation (12.621) et en calculant l'integrale angulaire, nous 
trouvons alors une equation de diffusion spatiale 

d t F(r;t)-^-V 2 r F(r;t) = (2.64) 

OA 

Le coefficient de diffusion (spatiale) est done donne par 

A P = ^- (2.65) 



29 



Discussion. Une situation physique interessante correspond a l'expansion d'un nu- 
age initialement confine et plonge dans un potentiel aleatoire. Apres un temps At, les 
atonies ont diffuse a travers une distance quadratique moyenne 

v 5 At 

Ar 2 ~ 6D sp At = — — ~ v 2 T iso At < (vAt) 2 (2.66) 
K 

A premiere vue, il semble que les atomes diffusent tres rapidement (dependance en 
t> 5 ). Mais la deuxieme formulation dans (12.661) nous rappelle que, lorsque le temps At 
est plus long que le temps caracteristique T iso de revolution angulaire, la diffusion est 
au contraire beaucoup plus lente que le vol libre. 

II reste a remarquer que pour un temps d'experience donne, ce sont les atomes 
les plus lents qui « ont le temps » d'entrer dans le regime diffusif (leur T iso oc v 3 est 
compris dans la duree de 1' experience). Les atomes plus rapides suivent par contre un 
transport balistique avec Ar 2 ~ (vAt) 2 , ils n'ont en outre pas encore « oublie » leur 
distribution angulaire initiale. C'est ainsi qu'ils ont parcouru la distance la plus grande 
(et que leur distribution en position permet de remonter a la distribution des vitesses 
initiale). Nous constatons done qu'il est possible d'observer simultanement plusieurs 
comportement du transport, pour chacune des classes de vitesses (radiales) des atomes. 
Le comportement diffusif apparait seulement pour les classes les plus lentes, il va se 
manifester, a la fin de l'experience, au centre de la distribution spatiale. 
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Chapitre 3 



Transport d'atomes — point de vue 
ondulatoire 

Un bref apercu d'approches ondulatoires pour modeliser le transport d'atomes 
dans un milieu aleatoire. Un grand nombre de resultats sont bien connus dans le 
contexte optique (« diffusion multiple de la lumiere»). Nous montrerons en partic- 
ulier comment Ton obtient une equation du transfert radiatif a partir de la descrip- 
tion ondulatoire. C'est un exemple particulier du passage du « microscopique » vers le 
« macroscopique ». 

3.1 Introduction 

Nous nous placons maintenant dans un contexte ondulatoire ou les atomes sont 
decrits, au niveau le plus fondamental de la theorie, par une fonction d'onde ip(r, t). 
Son equation devolution est l'equation de Schrodinger : 

ihd^(r,t) + ^V 2 ^ = V(r)^ (3.1) 

ou V(r) est le potentiel aleatoire. Nous nous interessons a des proprietes moyennes des 
atomes lorsque Ton ignore presque tout du potentiel aleatoire, sauf sa valeur moyenne 
(que nous supposons egale a zero, avec une redefinition convenable de l'energie 1 ) et 
sa fonction de correlation (V(ri) V(r 2 )). Les quantites interessantes dans le contexte 
de cette approche statistique sont, d'une part, la valeur moyenne (ip(r)) de la fonction 
d'onde et, d'autre part, safonction de correlation (ip(ri) ■0*(r 2 )) que nous appelerons 
egalement la fonction de coherence. Notons que la premiere de ces quantites n'est 

'Pour des atomes dans un potentiel aleatoire « profond » (qui peut devenir plus grand que l'energie 
cinetique atomique), cette redefinition de l'energie implique que la fonction d'onde moyenne est une 
onde « localisee » : son vecteur d'onde fc est imaginaire parce que (V) > Ei. Cette situation n'a pas 
d' analogue lumineux, et il serait interessant d'en explorer les consequences d'un point de vue theorique. 
Mais Ton n'aura pas le temps de le faire ici... 
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pas forcement nulle lorsque Ton se trouve au voisinage d'une « source coherente » 
d'atomes. La fonction d'onde moyenne decrit alors la partie « coherente » de l'onde 
atomique dont la phase est suffisamment bien definie pour qu'elle puisse interferer 
avec une onde de reference. La fonction de coherence, quant a elle, decrit la distribu- 
tion en position et en impulsion de l'ensemble atomique, comme nous verrons dans la 
representation de Wigner. Elle decrit egalement le contraste des franges lorsque Ton 
fait interferer les ondes rayonnees par deux points « du meme ensemble » (en opposi- 
tion a 1' interference avec une onde de reference dont il etait question pour la fonction 
d'onde moyenne). 

Dans la suite de cette introduction, nous allons fixer quelques notations. Quant au 
potentiel aleatoire, nous introduisons sa densite spectrale (spatiale) 

5(q) = Jdr (V(r + r„) V(r )) exp (-iq-r) (3.2) 

qui est la transformee de Fourier de la fonction de correlation du potentiel. Dans 
cette definition, nous supposons que le potentiel aleatoire est (statistiquement) ho- 
mogene, c'est-a-dire que sa fonction de correlation en deux points ri 2 ne depend que 
de la difference r 2 — ri. La portee de la densite spectrale (13.21) en fonction du vecteur 
d'onde q est de l'ordre de l'inverse l/£ c de la longueur de correlation. 
Quant a l'ensemble atomique, introduisons sa distribution de Wigner : 

/(r, k) = (V»(r + |s) V*(r - |s)) exp (-ik ■ s) (3.3) 

Cette distribution a des proprietes similaires a une densite de probabilite dans l'espace 
des phases : 

- elle depend a la fois de la position r et du vecteur d'onde k ; 

- son integrate sur les positions (sur les vecteurs d'onde) donne la distribution en 
vecteur d'onde (en position) des atomes ; 

- le flux moyen atomique peut s'exprimer comme la valeur moyenne du vecteur 
d'onde : 

Ti f /zk 

j(r) = — Im (V*(r) W(r)> = J dk— f(v, k); (3.4) 

- notons cependant que la distribution de Wigner n'est pas necessairement pos- 
itive, elle ne peut done pas en toutes circonstances etre interpreted comme une 
densite de probabilite. 

Nous allons souvent considerer des atomes avec une energie E fixee. L equation 
de Schrodinger stationnaire peut alors s'ecrire de la fa§on suivante : 

W(r)+A; 2 [l+/i(r)]V(r) = 0; (3.5) 
E = h 2 k 2 j2M, (i(r) = -V(t)/E. (3.6) 

Cette forme de l'equation met en evidence l'analogie a la lumiere (dans l'approxima- 
tion scalaire), ou /x(r) est la deviation de la constante dielectrique de l'unite. L'equation 
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de Schrodinger stationnaire sera notre point de depart pour la theorie de la « diffusion 
multiple des ondes de matiere». Une quantite importante dans ce contexte sera la 
fonction de correlation de la constante dielectrique, done du potentiel aleatoire renor- 
malise : 

( f i(r 1 ) f i(r 2 )) = E- 2 (V(r 1 )V(r 2 )) = e 2 g(r 2 - ri ), (3.7) 

e = (V 2 ) 1/2 /E (3.8) 

Le nombre e donne done l'ordre de grandeur du potentiel aleatoire par rapport a l'en- 
ergie des atonies. Au chapitre precedent, nous avons seulement considere le cas e< 1, 
Nous ferons de meme ici, notre theorie ne s' applique done pas au cas d'un potentiel 
aleatoire qui « piegerait » les atonies dans ses vallees. 

3.2 Point de vue de Ryzhik et ah 

Dans 1' article de Ryzhik, Papanicolaou et Keller [Keller96|, la diffusion 
des ondes scalaires par un potentiel aleatoire est etudiee en guise d' introduction. Nous 
nous contenterons ici d'esquisser leur demarche. 

3.2.1 Equation de transport a grande echelle 

lis cherchent une equation d' evolution pour la distribution de Wigner a des 
echelles spatiales grandes par rapport a la portee des correlations du potentiel aleatoire, 
d'une part, et a la longueur d'onde atomique, d' autre part. En utilisant un developpe- 
ment « multi-echelles », ils parviennent a demontrer que la distribution de Wigner 
verifie une equation du type transfert radiatif qu'ils ecrivent de la facon suivante : 

(* + £.V,)/(r, k ) = ^/^ S(k '- k )^- k V 

x[/(r,k')-/(r,k)] (3.9) 

Le dernier terme dans cette equation decrit la diffusion d'une onde atomique du vecteur 
d'onde initial k vers le vecteur d'onde final k'. L'efficacite de ce processus est pro- 
portionnelle a la densite spectrale du potentiel S(k' — k) au transfert de vecteur 
d'onde, comme e'est le cas dans 1' approximation de Born. L' equation du transfert 
radiatif (13.91) implique done les valeurs suivantes pour la section efficace de diffusion 2 
cr tot et la fonction de phase 3 p(n', n; ko) 

*tot(*b) = ^3 -—S(k (n'-n)) (3.10) 

Tin J 47T 

2 I1 s'agit en fait du nombre de diffusions par unite de temps. Elle est reliee au libre parcours moyen 
4x par crtot = hk / (M£ cx ). 

3 La fonction de phase (13 . 1 1 i depend des deux vecteurs unitaires n', n ; elle est normalisee telle que 
Jd 2 np(n',n;fc ) = 1. 
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p(n',n;k ) = k ° f S(k (n' - n)) (3.11) 

3.2.2 Discussion 

Tout d'abord, l'equation de transport (13.91 ) est un resultat utile pour fixer le com- 
portement a grande distance de la distribution de WlGNER. D'une facon plus con- 
ceptuelle, le travail de Ryzhik et al. fournit done une demonstration alternative de 
l'equation de transfert radiatif pour les ondes scalaires : celle-ci est valable a des 
echelles spatiales grandes par rapport a la longueur d'onde et la longueur de correla- 
tion. On en a vu d'autres demonstrations par Rytov ([Rytov|, § 4.3) et Baraba- 
NENKOV et Finkel'berg [Barabanenkov67 1. 

Ensuite, il ne faut pas oublier que (13.91) contient implicitement l'hypothese que 
le potentiel aleatoire est « faible » ; en effet, nous avons constate que la diffusion par 
le potentiel est decrite dans 1' approximation de Born (la section efficace est propor- 
tionnelle au carre de la transformee de Fourier du potentiel). Nous retrouvons la la 
proposition de LUCK HLuck93llLuck9 6l qui dit que l'equation de Bethe-Salpeter 
est equivalente a l'equation du transfert radiatif lorsque Ton tient compte des dia- 
grammes « en echelle » {ladder approximation) ce qui revient a une approximation de 
diffusion simple et done de Born. 

Finalement, nous remarquons que l'equation de transfert (13.91 ) se simplifie lorsque 
Ton se place dans le regime semi-classique. En utilisant la fonction de phase (13.1 IL 
nous constatons que la valeur maximale de Tangle de diffusion 9 est donnee par 

fc |n'- n| max = 2A; o sin(0 max /2) ~ j =^ 2sin(0 max /2) ~ ^ < 1 (3.12) 

Dans le regime semi-classique, la diffusion des atomes se produit done de preference 
vers l'avant, et nous pouvons utiliser la formulation de Fokker-Planck de l'equa- 
tion de transfert radiatif introduite au chapitre[T] Nous avons verifie que Ton trouve une 
equation de Fokker-Planck identique a celle du chapitreEJou nous nous sommes 
place d'emblee dans une description classique. 4 

3.2.3 Conclusion 

A des echelles spatiales grandes par rapport a la longueur de correlation du po- 
tentiel aleatoire et pour un potentiel aleatoire faible, la distribution de Wigner des 
atomes verifie une equation de transport du type transfert radiatif. La fonction de 
phase dans cette theorie s'exprime en fonction de la densite spectrale du potentiel 
aleatoire, comme e'est le cas dans 1' approximation de Born (de diffusion simple). 
Nous retrouvons exactement l'image classique du transport lorsque la longueur d'onde 

4 I1 faut a cet effet exprimer le coefficient de diffusion par la fonction de correlation du potentiel 
aleatoire, et utiliser un developpement a Fordre le plus bas en Xds/^c pour retrouver la fonction de 
correlation de la force aleatoire. 
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atomique est petite par rapport a la longueur de correlation du potentiel ; la distribution 
de WlGNER verifie alors une equation de Fokker-Planck. 

3.3 Diffusion multiple des ondes 

Introduisons maintenant un formalisme ondulatoire plus general qui est egalement 
en mesure de decrire la diffusion multiple des ondes atomiques. Nous allons presenter 
les objets et les equations de base d'une telle theorie, a savoir 

- la fonction de GREEN moyenne, qui permet de calculer la fonction d'onde 
moyenne et qui verifie l'equation de Dyson, 

- la fonction de coherence atomique et son equation d'evolution, l'equation de 
Bethe-Salpeter (B.-S.). 

Au paragraphe 13.41 suivant. nous etudierons de facon generale le lien entre l'equation 
de B.-S. et la theorie du transfert radiatif. Nous nous servons a cet effet de la represen- 
tation de WlGNER. 

Dans ce paragraphe, nous nous placons dans une situation stationnaire, l'energie 
des atomes est done fixee par le vecteur d'onde atomique k (valeur dans le potentiel 
moyen). 

3.3.1 Fonction de Green moyenne 

La quantite centrale pour decrire la fonction d'onde atomique est la fonction de 
Green G(r, r') pour l'equation de SchrOdinger stationnaire (13.51 ) : 

V*G(r, r') + kl [1 + fi(r)] G(v, r') = 6(r - r') (3.13) 

Cet objet decrit la fonction d'onde rayonnee par une source ponctuelle (de frequence 
E/h) situee a la position r' dans le milieu aleatoire. Pour une source spatialement 
etendue, nous trouvons la fonction d'onde en sommant les champs rayonnes par tous 
les points source. 

Equation de Dyson 

La fonction de Green (13.131 ) depend du potentiel aleatoire en chaque point, et 
il est impossible dans la pratique de la calculer explicitement. L'objet interessant est 
done la fonction de Green moyenne, la moyenne etant prise sur les configurations du 
potentiel aleatoire. Nous utiliserons la notation 

G{r - r') = (G(r, r')) (3.14) 

pour cette fonction de Green. Elle decrit done la valeur moyenne de la fonction 
d'onde rayonnee par une source ponctuelle. En ecrivant (13.141) . nous avons suppose 
que le milieu aleatoire est statistiquement homogene : la fonction de Green moyenne 
ne depend alors que la distance r — r' entre le point d' observation et le point source. 
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II s'agit maintenant de trouver une equation fermee pour la fonction de Green 
moyenne. II est utile a cet effet de transformer l'equation de Schrodinger (13.131) en 
une equation integrate : 

G(r, r') = G (r - r') + k 2 G J &v x G Q {r - n)//^)^, r'), (3.15) 

ou nous avons utilise la fonction de Green Go(r — r') pour l'espace libre : 

expi^olr-r'l 

G (r-r) = — — (3.16) 



L' approximation de Born consiste a resoudre (13.151) en remplacant au membre droit, 
G par la solution en espace libre G . On constate que cette procedure est le premier 
terme d'une solution par iteration. En allant au-dela de 1' approximation de Born, une 
telle solution genere des produits avec de plus en plus de facteurs /i(r). L'on peut alors 
prendre la valeur moyenne de cette serie, pour faire apparaitre les fonctions de cor- 
relation du potentiel aleatoire. En toute generalite, Ton aura affaire a des correlations 
d'ordre arbitrairement eleve, comme par exemple 

(/i(ri)/i(r 2 )//(r 3 )) 

Pour simplifier la theorie, nous allons supposer que le potentiel aleatoire a une statis- 
tique gaussienne; c'est-a-dire que les fonctions de correlations d'ordre superieur a 
deux peuvent s'exprimer au moyen de la fonction de correlation a deux points. Toutes 
les correlations a un nombre de points impair s'annulent alors, et pour les correlation 
a quatre points, par exemple, nous avons 

(/i(ri)/i(r 2 )//(r 3 )/i(r 4 ) = (/x(ri)/x(r 2 ))(/x(r 3 )/x(r4) + 

+ (M r i)M r 3))(M r 2)M r 4> + (^(ri)/i(r 4 ))(/i(r 2 )/i(r 3 ). 

La solution iterative reste quand meme complexe et contient une infinite de termes. 
Pour les organiser, Ton peut se servir d'une methode de diagrammes. Cette approche 
est exposee dans les articles de Frisch [Frisch66| et dans les livres de Rytov et al. 
([Rytov|, § 4.1) et de Ping Sheng ([Sheng95|, §§ 4.3, 4.4). En re-organisant la serie 



iterative, l'on trouve l'equation suivante pour la fonction de Green moyenne, l'equa- 
tion de Dyson : 

G(r - r') = G (r - r') + / d ri dr 2 G (r - r 1 )m(r 1 - r 2 )G(r 2 - r'), (3.17) 

ou m(ri — r 2 ) est appele l'« operateur de masse ». Son developpement iteratif ne con- 
tient que des « diagrammes irreductibles » 5 et s'ecrit : 

m(n - r 2 ) = eXg(ri ~ r 2 )G (r 1 - r 2 ) + (3.18) 

+ e 4 fc 8 |dr 3 dr 4 G'o(r 1 -r 3 )Go(r 3 -r4)G'o(r 4 -r 2 ) x (3.19) 

x [g(n - r 4 )5-(r 2 - r 3 ) + g(r 1 - r 2 )g(r 4 - r 3 )] + . . . 
5 " strongly connected" 
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Nous notons que l'operateur de masse ne depend que de la difference des positions 
i - ! — r 2 , ceci etant du a l'homogeneite statistique. 

Afin d'interpreter l'operateur de masse, revenons a une formulation differentielle 
de l'equation de Dyson, en appliquant l'operateur V 2 . + k\ a (13.171) : 

(v 2 + ko) G(r - r) = 5(r - r) + / dnm(r - n)^^ - r'), (3.20) 

Ceci est en fait une equation integro-differentielle qui exprime le champ moyen rayon- 
nee dans le milieu aleatoire comme la somme d'un champ en espace libre (la fonction 
S), plus une correction qui depend de facon non-locale de la valeur du champ autour 
de la source (le deuxieme terme). L'operateur de masse traduit done precisement la 
retro-action du milieu diffusant sur la propagation du champ moyen. 

Nous notons finalement que Ton peut resoudre (13.201 ) par une transformee de 
Fourier puisque le milieu est homogene. L'integrale du deuxieme terme est en effet 
un produit de convolution, et nous avons done pour les transformees de Fourier 

(£; 2 -k 2 )G(k) = l + m(k)G(k) (3.21) 

d '°- : 5 < k > = w^hm (3 - 22) 

Nous en deduisons la relation de dispersion des ondes de matiere dans le milieu desor- 
donnee : 

k 2 + m(k) = k 2 (3.23) 

Cette equation definit le vecteur d'onde « effectif » k C R de la fonction d'onde moyenne 
dans le milieu. II differe du vecteur d'onde dans le vide k pour deux raisons : d'une 
part, la diffusion par le milieu introduit de 1' attenuation pour le champ coherent, ce 
qui se traduit par une partie imaginaire du vecteur d'onde, et d' autre part, la diffusion 
modifie la phase de propagation des ondes. 6 



3.3.2 Approximation de Bourret 

Tout ce que nous venons d'exposer ne serait en fait que des manipulations formelles 
d'equations s'il n'y avait pas un moyen de calculer l'operateur de masse m(r 1 — r 2 ). 
Or, le developpement (13.191 ) permet de le faire en principe (il contient neanmoins une 
infinite de termes), et en particulier dans une limite perturbative, lorsque Ton suppose 
f«l.A l'ordre le plus bas, Ton trouve done l'expression assez maniable de la pre- 
miere ligne (13.181) 

m(n - r 2 ) w e 2 k 4 g{v x - r 2 )G (r 1 - r 2 ), (3.24) 

6 Rytov dit que la longueur des rayons augmente par l'effet de la diffusion et par consequent aussi 
la partie reelle de k e ff ( |Rytov| , § 4.2, p. 137). 
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que Ton appelle, suivant les auteurs, l'« approximation de BOURRET » HFrisch66L 
l'« approximation de regularisation du premier ordre » [Frisch66|, la "ladder approxi- 
mation" [ |Rytov| . Elle doit egalement avoir un lien avec la "random phase approxima- 
tion" (RPA) utilisee dans la physique des solides. 

Nous remarquons que 1' approximation de BOURRET (13.241) ressemble, au premier 
abord, a 1' approximation de Born pour l'operateur de masse. Qu'a-t-on alors gagne 
par rapport au choix de faire d'emblee 1' approximation de Born pour la fonction de 
Green ? La reponse se trouve par exemple dans la formule (13.231) qui permet de cal- 
culer les « nouveaux modes » dans le milieu aleatoire, alors qu'avec 1' approximation 
de Born, la fonction de Green donne seulement une amplitude de diffusion. Da 
facon plus profonde, le formalisme de l'equation de Dyson permet, dans 1' approxi- 
mation de Bourret, de trouver une fonction de Green moyenne qui, elle, est deja 
une resommation partielle d'une infinite de diagrammes. 

Pour illustrer la puissance de 1' approximation de BOURRET pour la diffusion mul- 
tiple des ondes, citons le resultat pour le vecteur d'onde effectif dans un milieu faible- 
ment diffuseur, donne par RYTOV (| |Rytov| |, eq.(4.61)) : 



^ = ^ + ^/^ 9 OTlog(^i|;) 2 + i^7^ 9 W (3.25, 

ou la fonction de correlation est supposee isotrope (sa transformee de FOURIER g(k') 
ne depend alors que du module kl du vecteur d'onde). On constate que le vecteur 
d'onde dans le milieu contient une partie imaginaire positive (attenuation de l'onde 
coherente par la diffusion), et une partie reelle un peu plus grande que dans le vide. 

Pour donner un ordre de grandeur de ces effets, prenons un potentiel aleatoire avec 
une correlation gaussienne. La densite spectrale est alors gaussienne egalement. La 
transformee de Fourier m(k) de l'operateur de masse est alors donnee par 

m ( k ) = ^(efcglc) 2 [ e -(fco+fc)^/2 erfc[ _ i(A;o + fc)4 /v /2] (3.26) 

2 f°° 9 

erfc(ar) = -= / e~ l dt (3.27) 

\/7T Jx 



Cette fonction est representee sur la figure l3~Tl Quant a la relation de dispersion (13.231) . 
Ton trouve le vecteur d'onde effectif suivant (dans la limite semi-classique k £ c 3> 1) 

hs « k (l + Yj + ifco 2(2 ^ 3/2 e 2 fc 4- (3.28) 

Rytov n'oublie pas de nous rappeler que ce resultat n'est valable que dans l'approxi- 
mation de BOURRET. En comparant aux corrections d'ordre suivant, Rytov trouve la 
condition de validite suivante 

BOURRET : (ek Q 2 < 1. (3.29) 
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FIG. 3.1 - Operateur de masse pour un potentiel aleatoire a correlations gaussienne 
isotropes. Trait plein : partie reelle ; tirets : partie imaginaire. 

m(k) est donne en unites de k$, la longueur de correlation vaut i c = 10/ ko, et le petit 
parametre de 1' approximation de BOURRET vaut (ek £c) 2 = 0.1. 



Rappelons que dans les methodes classiques du chapitre El la condition de validite 
d'une approche de Fokker-Planck etait que le potentiel soit faible par rapport a 
l'energie cinetique des atomes, ce qui revient dans la notation de ce chapitre a la con- 
dition e 2 < 1 [voir (12.331) et (I2.35l) l. L' approximation de BOURRET (13.291 ) impose 
done une limite plus stricte au potentiel moyen qui Test d'autant plus que Ton se place 
dans le regime semi-classique ou k £ c 1. 



Remarque. L' approximation de BOURRET peut etre generalisee a un mi- 
lieu aleatoire dont la statistique n'est pas gaussienne. Les fonctions de cor- 
relations peuvent alors etre ecrites comme une somme sur des correlations 
« elementaires » gk(?i, ■ ■ ■ , ffe) pour un groupe de k points (developpement en es- 
saims 7 ). Finkel'berg montre alors que pour calculer l'operateur de masse a l'or- 
dre le plus bas, il suffit de retenir les diagrammes irreductibles qui ne contiennent 
qu'un seul groupe de points |Fihkelberg67]. II appelle cette approximation la "single- 
group approximation". Son avantage est qu'elle permet egalement de tenir compte de 
la diffusion par des agglomerats de plus que deux particules. Finkel'berg donne la 
condition de validite suivante pour cette approximation : 



single group : 



dm 



dk 2 



< 1. 



(3.30) 



''"cluster expansion'' 
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II montre en particulier que cette condition est suffisante pour que la diffusion soit 
faible (libre parcours moyen beaucoup plus grand que la longueur d'onde). 

3.3.3 Equation de Bethe-Salpeter 

Pour finir ce paragraphe, il ne reste qu'a presenter l'equation d'evolution pour la 
fonction de correlation de la fonction d'onde atomique. Nous allons a cet effet con- 
siderer non plus la fonction de Green, mais l'onde rayonnee par une source g(r) 
quelconque 

V>(r) = fdr'G(r,r r )e(r') (3.31) 

L'onde moyenne (ip(r)) s'exprime alors par une equation analogue ou intervient la 
fonction de Green moyenne G(r — r'). En ce qui concerne la fonction de coherence 
des atonies, elle s'exprime par une correlation entre deux fonctions de Green : 

MrOV^ra)) = J dr[ dr 2 (G(r u r[)G*(r 2 , r' 2 )> 6 (r[) 6 * (v' 2 ) (3.32) 

Cette fonction de correlation entre deux fonctions de Green est l'objet de l'equation 
de Bethe-Salpeter. Elle se deduit de l'equation de Schrodinger (13.51) de facon 
analogue a l'equation de Dyson (13.171) (voir HFrisch 66 1 et § 4.3 de [Sheng95| pour 
plus de details). L'equation de Bethe-Salpeter permet d'ecrire l'equation fermee 
suivante pour la fonction de coherence du champ atomique : 

(V(ri)V*(r 2 )) = ^(r 1 ))(^(r 2 )) + |dr / 1 dr / 2 dp 1 dp 2 (3.33) 

G( ri - r;)cT(r 2 - r 2 ) K(v[, Pl \v' 2 , p 2 ) (^(Px)^(p 2 )) 

Dans cette equation apparait la fonction K(r 1 , p 1 \r 2 , p 2 ) que Frisch appelle 
l'« operateur d'intensite ». Chez d'autres auteurs, elle porte le nom de « vertex irre- 
ductible » [Tsang85 , Sheng95 1. D'un point de vue physique, cet objet decrit la correla- 
tion aux positions r 1]2 entre les champs rayonnes par deux sources situees a p 12 . Son 
developpement perturbatif (en diagrammes) contient les premiers termes suivants : 

^(ri,Pi|r 2 ,P 2 ) = e 2 k^g(r 1 - r 2 )5(r 1 - Pl )8(r 2 - p 2 ) + (3.34) 
+ e% g{r x - r 2 ) g(p 1 - p 2 ) [5(n - Pl ) C7 (r 2 - p 2 ) + {1 <-> 2}] + 
+ e 4 ^o 9iji ~ P 2 ) g{pi - r2)G (ri - Pl ) C7 (r 2 - p 2 ) + 
+ eX g( ri - Pl ) 5(r 2 - p 2 ) J dr\ G ( ri - r[) g(r[ - r 2 ) G (r' 1 - Pl ) + 
+ {1^2} + ... (3.35) 

L'operateur d'intensite est done un objet plutot encombrant. Son approximation de 
BOURRET est donnee par la premiere ligne, l'eq. (13.341) . ou Ton constate qu'il revient 
essentiellement a la fonction de correlation du potentiel si ce dernier est faible. 
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Notons encore la propriete generale suivante de l'operateur (13.351) . elle aussi due a 
l'homogeneite du milieu diffuseur : il est invariant par une translation globale de ses 
quatre arguments 

K(t! + Ar, Pl + Ar|r 2 + Ar, p 2 + Ar) = K(t u Pl \r 2 , p 2 ) (3.36) 

Par consequent, il ne peut dependre que des differences des positions, comme on le 
constate aussi dans le developpement (I3.35I) . Ceci se traduit par la propriete suivante 
pour la transformee de Fourier de l'operateur d'intensite : 

K(k lt «i|k 2 , k 2 ) = J dri dr 2 dp 1 dp 2 

K(t u p 1 |r 2 ,p 2 )expi (-ki-ri + K V p 1 + k 2 r 2 - K 2 -p 2 ) (3.37) 
= (2tt) 3 <J(ki - ki - (k 2 - k 2 )) K(ki, Kilka, k 2 ) (3.38) 

Les signes des vecteurs d'onde ont ete choisis tels que la T.-F. de l'operateur d'intensite 
admet 1' interpretation suivante : il decrit la correlation entre les amplitudes pour les 
deux processus de diffusion 

Ki — > ki et k 2 — ► k 2 

La fonction 5 dans (13.381) exprime alors que ces deux processus ne sont correles que 
lorsqu'ils font intervenir le meme transfert de vecteur d'onde : 

K ^ <J=^ ki - ki = k 2 - k 2 (3.39) 

Cette propriete se comprend aisement dans 1' approximation de Born : 1' amplitude de 
diffusion pour le processus k, x — > ki est proportionnelle a la composante de Fourier 
/i(ki — ki) du potentiel aleatoire, et les composantes de Fourier du potentiel aleatoire 
sont « 5-correlees » lorsque le potentiel est statistiquement homogene : 

(MqiKM) = S(qx - q 2 )(27r) 3 e 2 ( 7 (q 1 ). (3.40) 

Dans la pratique, il n'est pas possible de resoudre l'equation de Bethe- 
Salpeter (13.331) . meme dans 1' approximation de BOURRET. Au paragraphe suivant, 
nous retra5ons une reformulation de l'eq. de B.-S. qui permet de retrouver l'approche 
du transfert radiatif. 

3.4 Lien a l'equation du transfert radiatif 

L'interet de ce paragraphe est de preciser le domaine de validite de l'equation du 
transfert radiatif, en la deduisant du formalisme microscopique de la diffusion mul- 
tiple des ondes. Une presentation semblable a ete donnee par Barabanenkov et 
Finkel'berg IBarabanenk ov67l et par RYTOV et al. ([Rytov|, § 4.3) pour des on- 
des scalaires. Pour des ondes electromagnetiques, une generalisation de l'equation du 
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transfert radiatif (portant sur le vecteur de Stokes) est deja formulee dans le livre de 
Chandrasekhar. Une justification micro scopique en a ete donnee par K. M. WAT- 
SON HWatson 69l dans une situation de diffusion faible, et J.-M. Tualle a fourni une 
presentation plus generale [TualleT|. Une bonne reference est egalement 1' article de 
revue de Lagendijk et van Tiggelen jvanTiggelen96). 

Nous verrons que 1' equation du transfert radiatif (ETR) est valable a grande echelle 
spatiale (ce qui n'est guere etonnant, vu le resultat de Ryzhik et ah), mais aussi 
quelle que soit la force du potentiel diffuseur. Ce dernier point semble ne pas avoir 
ete souleve par d'autres auteurs parce que Ton se place dans la plupart des cas dans le 
regime de diffusion faible, en utilisant 1' approximation de BOURRET. Luck, par ex- 
emple, dit dans ses papiers HLuck93[ ILuck96l . au tournant d'une phrase : « l'equation 
de Bethe-Salpeter, dans 1' approximation de Bourret, est equivalente a l'equa- 
tion du transfert radiatif ». Nous allons par contre rester assez general et ne pas faire 
1' approximation de BOURRET. Ainsi le formalisme permet-il de voir a quel endroit 
l'equation du transfert radiatif est moins generale que celle de B.-S. 

Nous allons seulement faire l'hypothese suivante (qui ne semble pas tres restric- 
tive) 

- le milieu est statistiquement homogene. La fonction de Green moyenne G, 
qui intervient dans l'equation de B.-S. (13.331 ). admet alors le developpement de 



ou m(k) est l'operateur de masse. 

3.4.1 Bethe-Salpeter dans Wigner 

La premiere etape consiste a exprimer l'equation de B.-S. (13.331) en representation 
de Wigner. En introduisant les transformees de Fourier des fonctions de Green 
moyennees et de l'operateur d'intensite, Ton trouve l'equation suivante : 



Calculons d'abord le premier terme qui correspond a la transformee de Wigner 
du produit des champs moyens. Nous exprimons le champ moyen a l'aide de la 
fonction de Green moyennee (cf. eq. 13.311) et inserons ensuite la decomposition de 
Fourier (13.411) de cette derniere : 



Fourier (l3~2"2l 



,8 




(3.41) 



/(r, k ) = /coh(r, k) + /as (r, k) 



(3.42) 



/coh(r,k) = 



8 Nous nous servons de la notation abregee 



dk = 



dk 
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= Jdse iks j & Pl &p 2 G(r + \s - Pl )G*{r-\s- p 2 ) g( Pl ) g*(p 2 ) 
= |dsdk 1 dk 2 dpG(k 1 )G*(k 2 )/ e [p,i(k 1 + k 2 )] x 

x exp i [-k- s + | (kj + k 2 ) • s + (k x - k 2 ) • (r - p)] (3.43) 

Nous avons effectue une des deux integrations sur les positions de la source g(p) pour 
faire apparaitre sa transformee de Wigner f e [p, |(ki + k 2 )]. L' integration sur s est 
maintenant immediate et l'equation (13.431) devient (q = ki — k 2 , k = |(ki + k 2 )) : 

/ coh (r,k) = ydpdqi(q;k)/ e (p,k) expiq-(r-p) (3.44) 

avec 1' abbreviation suivante pour le produit des transformees de Fourier des fonc- 
tions de Green moyennees : 

i(q;k) = G(k+±q)G*(k-±q) (3.45) 

Nous constatons done que le terme coherent relie la luminance observee a celle de la 
source au meme vecteur d'onde k. Ecrivons ce resultat encore sous la forme integrale 
d'une equation de transfert radiatif. Le terme coherent s'exprime comme une integrale 
sur les points sources, ponderes avec une « fonction d' attenuation » (la transformee de 
sc Fourier de A(q; k)) : 

/ coh (r,k) = JdpA(r-p;k)f e (p,k), (3.46) 
A(r - p; k) = J dqi(q; k) exp iq- (r - p) (3.47) 

Pour le deuxieme terme f^m dans ( 13.42b . nous obtenons d'abord, par une procedure 
similaire, l'expression suivante qui ressemble a (13.431) : 

/diff(r, k) = J ds dki dk 2 <\p&n x &k 2 

K(k u Kl |k 2 , k 2 ) f[p, + k 2 )] G(k x ) G*(k 2 ) x 

xexpi -k-s+ |(ki + k 2 )-s+ (ki - k 2 )-r - (^i - K 2 )-p (3.48) 

Les integrations sur s et sur la demi-somme k = ~(ki + k 2 ) des vecteurs d'onde s'ef- 
fectuent de la meme fa§on qu'en (13.431) pour faire apparaitre la fonction A(q; k) (I3.45I) . 
Pour l'integration sur «i )2 , passons a la mi-somme et la difference, et utilisons l'invari- 
ance par translation (13.381) de l'operateur d'intensite. Nous trouvons ainsi l'expression 
assez compacte similaire a (13.461) (avec k = + k 2 )) : 

/difl(r, k) = J dp$Kp A (r - p; k, n)f(p, k) (3.49) 
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oil intervient une fonction d' attenuation differente qui ressemble plutot a une fonction 
de phase (avec q = k x — k 2 = n\ — n 2 ) '■ 

p A (r- p; k, k) = 

J dq exp[iq- (r - p)] I(q; k) k(k + |q, k + |q|k - |q, K - |q) (3.50) 

Ce qui est remarquable, c'est que nous n'avons encore fait aucune approximation pour 
deduire ce resultat. 

3.4.2 Comparaison a l'equation de transfert 

Regardons d'abord ce que Ton a pu faire : l'equation de Bethe-Salpeter a ete 
transformee en : 

/(r, k) = J dpA(r - p, k)/ ff (p, k) + JdpSKp A (r - p; k, rc)/(p, k). (3.51) 

La distribution de Wigner (rappelons qu'elle est equivalente a la luminance pour 
la lumiere) contient un premier terme « coherent » (13.461 )) qui correspond aux champ 
moyen rayonne par les sources dans le milieu. Le deuxieme terme, l'integrale (13.491) . 
doit alors correspondre aux ondes diffusees par le milieu. II a bien la structure d'une 
integrale de diffusion, ou intervient la fonction de Wigner pour un vecteur d'onde k, 
different. En effet, le vecteur d'onde k peut tres bien s'interpreter comme un « vecteur 
d'onde entrant » qui est diffuse vers le vecteur d'onde k, avec une probabilite donnee 
par l'operateur d'intensite. 

Les integrales sur p dans les deux termes de (13.511) ne sont pas si surprenantes 
qu'elles le semblent au premier abord, parce que nous sommes en train d'etudier une 
formulation integrale, qu'il faut comparer a la version integrale de l'equation du trans- 
fert radiatif. En utilisant la forme donnee par Chandrasekhar ( [ Chandrase kharl . 
Chap. 1, eq.(50)) et en ajoutant la luminance S(r, n) des sources, la version integrale 
de l'ETR s'ecrit (/ est la luminance) : 

oo oo 

J(r, n) = e" s/ ^(r -sn,n)+ J y e" s ^ J dn'p(n , ) n) J(r - sn, n') (3.52) 

o o 

ou r — sn est un point sur le rayon de direction n a distance s du point d'observa- 
tion. L'integrale sur s exprime le fait que la luminance en r est le produit de tous les 
evenements de diffusion survenus aux positions anterieures a r sur ce rayon et qui ont 
amene la lumiere dans la direction n. Par ailleurs, le premier terme dans (13.521) a bien 
la forme de la partie « coherente » du champ, qui est attenuee par la diffusion sur une 
distance caracteristique egale au libre parcours moyen i. 

A la difference du transfert radiatif (I3.52I ). l'integrale sur p dans l'expression ex- 
acte (13.511 ) s'etend a priori sur tout l'espace. La fonction de Wigner verifie done une 
equation de transport non-locale. Pour retrouver l'ETR standard (13.521) qui, elle, est 
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locale, c'est l'integrale sur q dans la fonction d' attenuation A(r — p; k) (13.471) qui doit 
limiter l'integrale sur p de sorte qu'elle ne porte plus que sur le rayon arrivant au point 
d'observation r de la direction k. Nous allons maintenant preciser les conditions pour 
qu'apparaisse une telle simplification. 



3.4.3 Solution approchee a grandes distances 

Astuce 

Transformons d'abord le produit des deux denominateurs dans (13.451) selon la for- 
mule suivante (evidente, mais astucieuse. . .) 

i(q;k)=G(k+iq)G*(k-iq) 



[kl - (k + |q)2 - m (k + |q)] [k$ - (k - |q)a - m*(k - ±q)] 

(3.53) 



G(k+^q)-G*(k-^q) 



2kq + m(k + |q) — m*(k — |q) 

Cette expression est encore exacte. Si nous cherchons maintenant le comportement de 
la distribution de WlGNER a des echelles |r — p\ grandes par rapport a la longueur 
d'onde Xjb, l'integrale (13.471) montre que Ton peut se contenter des vecteurs d'onde q 
tres petits par rapport aux vecteurs d'onde atomiques k, k. L'idee vient alors a l'esprit 
de faire un developpement limite en fonction de q. Nous allons le faire ici a l'ordre le 
plus bas ; Ton verra que Ton retrouve alors la forme (13.521) de l'equation du transfer! 
radiatif. Au paragraphe 13.4.41 suivant, nous irons un ordre plus loin pour obtenir une 
premiere correction par rapport a l'equation de transfert. 
A l'ordre le plus bas en q, la formule (13.531) devient 

i (q;k ). ^g-^ , (3,4, 
2k- q - 2ilm/4 r 

ou nous avons utilise le vecteur d'onde effectif (complexe) pour ecrire le pole du 
denominateur : 

m(k ± |q) « m(k) « jfcg - k 2 eS (3.55) 

Nous avons egalement approximee la partie imaginaire de la fonction de Green 
moyennee (le terme au numerateur dans (13.531) ) par une fonction 5. Celle-ci exprime 
le fait que le module des vecteurs d'onde atomiques finaux k (apres le processus de 
diffusion) est fixe par la relation de dispersion dans le milieu (la partie reelle du vecteur 
d'onde effectif k e g). Cette approximation est justifiee dans le regime de faible diffusion 
ou le libre parcours moyen est beaucoup plus grand que la longueur d'onde. Une telle 
situation correspond a l'image intuitive que Ton se fait de l'ETR. La limite opposee 
correspond a des ondes localisees dans le milieu aleatoire, un regime qui demande une 
description au-dela de l'ETR standard. 
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Retrouver l'ETR differentielle 



En reportant (13.541) dans l'ETR et en negligeant le petit vecteur d'onde q dans 
l'operateur d'intensite, nous observons que l'ETR (13.511) se resoud facilement par une 
transformee de Fourier par rapport a la position r. Nous notons /(q, k) la fonction 
de Wigner transformee, et elle verifie 1' equation suivante 

(ik-q + Im/4) /(q,k) = 
7r<y(fc 2 -ReA&) (/ fl (q,k) + J dnK(k, *|k, K )/(q, «)) (3.56) 

En revenant dans l'espace directe, nous trouvons done la forme differentielle de l'ETR, 
dont le libre parcours moyen et la fonction de phase sont donnes par 



1 I m ^eff Imm(k p 



(3.57) 



p(u,u') = ( 4?rfc y K(kp u i k p u'\k p u, k p u') (3.58) 

k p = (Re^) 1 / 2 (3.59) 

Nous avons done trouve ici des expressions microscopiques pour les parametres qui in- 
terviennent dans l'ETR. Des formules identiques ont etes ecrites par Barabanenkov 
et FlNKEL'BERG (eq. (22) de IBarabanenk ov67l '). mais en utilisant une forme ap- 
prochee pour les operateurs de masse et d'intensite (approximation de Bourret). 
La demarche presentee ici montre que ces expressions restent valables dans un con- 
texte plus general, a condition que les ondes ne soient pas localisees dans le milieu. La 
seule difference par rapport aux resultats de BARABANENKOV et FlNKEL'BERG est 
1' apparition du vecteur d'onde k p qui decrit l'indice de refraction du milieu aleatoire. 
A la limite perturbative, il faut en effet prendre k p = k pour etre consistent avec la 
conservation de l'energie (l'identite de Ward, voir [vanTiggelen96|). 



Retrouver l'ETR integrale 

Nous pouvons egalement retrouver la formulation integrale de l'ETR donnee 
en (I3.52I ). A cet effet, nous calculons l'integrale sur q dans (13.511 ) dans 1' approxima- 
tion (13.541 ) pour la fonction d'attenuation A(q; k). L'integrale a d'ailleurs ete calculee 
au chapitre precedent (voir la note|3]en bas de la page^U) : 

exp iq-(r-p) 

dq 2k.q-2ilm^ ff = (3 - 60) 

(i/2k)8(v ± - p ± ) exp[-(Im k 2 eS /k 2 )k-(r - p)], 

lorsque k- (r — p) > 0; 
0, lorsque k- (r — p) < 0. 
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(L'indice « _L » fait reference au vecteur k.) Nous en concluons que l'integrale sur 
p ne porte en fait que sur un demi-rayon qui aboutit au point d' observation r en se 
propageant le long de la direction k. Le long de ce rayon, l'intensite du champ est 
attenuee avec une longueur caracteristique 

un fi, eff 

ce qui correspond bien a 1' attenuation du champ coherent par diffusion. Le terme co- 
herent s'ecrit done comme une integrate sur les sources « interceptees » le long de ce 
rayon, qui rayonnent dans la direction k : 

oo 

/ coh (r, k) = 5(k 2 - Re k 2 cS ) J ds' e~ s> ' 1 f e (r - s'k, k) (3.62) 

o 

Nous retrouvons ici le premier terme de l'equation du transfert radiatif de Chan- 
drasekhar (13.521) . La fonction 5 assure que les seuls vecteurs d'onde de la source 
qui rayonnent sont ceux qui correspondent a la relation de dispersion dans le milieu 
aleatoire. (Rappelons encore que la source g(r) dans (13.311) et (13.621) n'a pas la meme 
dimension que le champ, ceci est a l'origine du prefacteur.) 

Quant au deuxieme terme de (13.521) . qui decrit la diffusion, nous remarquons qu'a 
l'ordre le plus bas, l'operateur d'intensite dans (13.501) ne depend pas du vecteur d'onde 

K(k + ~q, k, + |q[k - |q, n - w K(k, «|k, k) (3.63) 

On peut alors le sortir des integrales sur q et p et definir une fonction de phase locale 
comme suit (attention, sa normalisation est differente de (13.581) ) 

p A (r - p;k, k) = A(r - p; k) p(k, k) 

p(k,«) = T$(k 2 -Rek* s )K(k,K\k,K). (3.64) 

k 

En utilisant le resultat (13.601) pour la fonction d' attenuation A(r — p; k), le terme de 
diffusion ( 13.491) prend done la forme 

oo 

/di ff (r,k) = Jds'e- s ' /e JdKp(k,K)f{r- 8 %K) (3.65) 
o 

Puisque la relations de dispersion fixe le module du vecteur d'onde k, le terme de 
diffusion est exactement celui donne par Chandrasekhar en (I3.52I) . (La fonction de 
phase (13.641) n'est pas normalisee de la meme facon que celle de Chandrasekhar, 
ceci explique la difference entre les prefacteurs.) 

Nous constatons done qu'aux grandes echelles spatiales, le transport de la distri- 
bution de Wigner des atomes (ou du champ electro-magnetique) est regi par une 
equation de transport radiatif. Ceci est vrai quelle que soit la force du potentiel alea- 
toire. Pour un potentiel faible, 1' approximation de BOURRET (13.341) de l'operateur 
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d'intensite montre que la fonction de phase est proportionnelle a la densite spectrale 
du potentiel aleatoire 

BOURRET : K(k + ~q, k + |q|k - ~q, k - |q) = 

AM 2 

= e 2 k* g(k -k) = —^S(k - K ) (3.66) 

n 

Nous retrouvons alors l'approche de Ryzhik et ah, ainsi que l'idee que s'est faite 
Luck de l'equation de transfer!. La demonstration que nous venons de donner de 
cette equation montre cependant qu'elle est egalement valable pour des potentiels plus 
forts ou 1' approximation de BOURRET n'est plus possible. 



Conclusion 

Nous avons constate qu'aux grandes echelles spatiales, la theorie de la diffusion 
multiple des ondes prend la forme d'une equation du transfert radiatif. Ce resultat 
repose seulement sur l'hypothese d'un milieu invariant par translation et de la diffusion 
faible (libre parcours moyen beaucoup plus grand que la longueur d'onde). Dans ce 
regime, la transformee de Wigner de la fonction de correlation decrit la luminance du 
champ, et elle verifie une ETR ( I3.56I ) dont les parametres macroscopique s'expriment 
a l'aide des operateurs de masse et d'intensite. 

Nous avons obtenu l'ETR ainsi generalisee par un developpement a l'ordre le plus 
bas par rapport au vecteur d'onde q, la variable conjuguee au profil spatial en r de la 
luminance. Cette procedure (la « limite de Kubo » [vanTiggelen96|) traduit formelle- 



ment la limite d'une « grande echelle spatiale». Au paragraphe suivant, nous allons 
pousser le developpement jusqu'au second ordre en q, afin de preciser l'echelle spa- 
tiale minimale au-dela de laquelle l'equation du transfert radiatif est justifiee. 



3.4.4 Echelle spatiale minimale 

pour l'equation du transfer radiatif 

Dans le calcul de la fonction de Wigner du champ aux grandes distances, nous 
avons pour l'instant seulement retenu les termes a l'ordre le plus bas en q. Nous al- 
lons ici aller plus loin et analyser les termes de l'ordre suivant pour trouver une limite 
superieure g max a ce developpement. Dans l'espace reel, la valeur g max se traduit par 
une echelle spatiale minimale L min = l/g max , avec la signification physique la suiv- 
ante : il faut que la fonction de Wigner du champ, en fonction de la position r, varie 
peu a l'echelle L mhl pour que l'equation du transfert radiatif soit valable. 

Nous allons nous contenter d'etudier la dependance de q de la fonction d' attenua- 
tion A(q; k) (13.471) . Quant a la variation de l'operateur d'intensite K ( 13.631 ) avec q, Ton 
en trouvera une discussion dans l'article de revue de KRAVTSOV et APRESYAN (§ 6.1 
de HKravtsov 96 1) : elle peut etre reliee a l'effet de memoire de la figure de tavelures 
lorsque Ton fait varier la direction d'incidence du faisceau. Nous notons egalement 
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que pour un potentiel aleatoire faible (approximation de Bourret), la formule (13.661) 
montre que K est independant de q. Un moyen d'apprecier sa dependance de q est 
alors de calculer des termes superieurs (13.351) qui vont au-dela de 1' approximation de 
Bourret. 

Poussons done le developpement de la fonction d'attenuation A(q; k) jusqu'a l'or- 
dre quadratique. Nous allons negliger la variation avec q de la difference des fonctions 
de Green 

G(k + |q) - G*(k - §q) « -2m5(k 2 - Re k 2 eS ) 

et seulement analyser le denominateur dans ( 13.531) qui presente un pole en q. A l'ordre 
quadratique, celui-ci devient : 

2k q + m(k + lq) - m*(k - §q) « 

2k- q - 2i Im k 2 s + q- Re V ' k m + - V q i q j Im (3.67) 

4 ^ dkidkj 

Les derivees de l'operateur de masse m(k) sont prises ici pour k verifiant la relation 
de dispersion. Pour continuer le calcul, nous allons faire rhypothese supplementaire 
suivante : 

- le milieu aleatoire est statistiquement isotrope. La relation de dispersion s'ecrit 
alors k = k e s et l'operateur de masse m(k) depend seulement du module k du 
vecteur d'onde k. 

Par exemple, le milieu a fonction de correlation gaussienne dont nous avons donne 
l'operateur de masse en (13.26b . est isotrope. Par consequent, les derivees de l'operateur 
de masse de simplifient : 

- dm 

V k m = k— 
ak 

d 2 m . 1 1 dm\ ? t ( d 2 m 1 dm\ 

r-rr l+hkj -rrr-T-rr (3-68) 



dhdkj 13 \k dk J 3 \ dk 2 kdk / 
Le denominateur (13.671) peut alors etre ecrit sous la forme suivante : 

2k-q + m(k + |q) -m*(k- |q) w 2k'-q - i 7 + iaq 2 - i/3(k-q) 2 (3.69) 

k '- k ( 1 + ^ Re 3 <3 ' 70) 

7 = -2lmm{k cS ) = 2lmk 2 s (3.71) 

La derivee premiere de l'operateur de masse change done simplement le module du 
vecteur d'onde k i— > k'. Le derivee seconde ajoute des termes quadratiques en q. 
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Avant de continuer le calcul, donnons une estimation jusqu' a quel vecteur d'onde 
maximal ce developpement est justifie. La figure l3~T1 et l'expression (13.261) montrent 
que l'operateur de masse m(k) varie sur une echelle caracteristique de l/£ c en fonction 
de k. Nous trouvons done la condition de validite 

|q| £ j (3-74) 

Dans l'espace reel, cette condition traduit le fait que dans le formalisme du transfert ra- 
diatif, la luminance ne resout pas la structure microscopique (longueur de correlation) 
du milieu diffusant. 



Calcul de la fonction d'attenuation 

Nous avons done 1' integrate suivante a calculer : 

exp iq-r 



/ 



dq 



2k' q — v) + ictq 2 — i/5(k-q^ 



2 



fTOrT expi(Q-R + gz) 

H q 2k' q - i 7 + iaQ 2 + i(a - Pk 2 )q 2 ' K J 



Nous avons decompose les vecteurs q et r en composantes (Q, q) , (R, z) perpendic- 
ulaires et paralleles au vecteur k'. La coordonnee z mesure done la distance du point 
d' observation d'un point source situe sur le rayon qui arrive de la direction k au point 
d' observation. 

II se trouve que Ton peut discuter de fa§on analytique le comportement de l'inte- 
grale (I3.75I) . Nous prenons d'abord l'integration sur la composante longitudinale q, en 
utilisant le theoreme des residus (nous avons abrege 7 = 7 — aQ 2 , a = a — (3k 2 ) : 

l &( i ou eXPi9 !- 2 = -r SgnZ s [©(^lmgi)e^ - 0(^lmg 2 )e^-] (3.76) 
J 2k'q — 17 + iaq z a\q\ — q 2 > 1 J 

ou Q(z) est la fonction de marche d'escalier et les deux quantites imaginaires gi 2 sont 
les solutions d'une equation quadratique 



(3.77) 



(Les signes sont choisis tels que < \q%\.) 

Pour interpreter et simplifier ce resultat, supposons que le potentiel aleatoire est 
faible (approximation de Bourret). L'operateur de masse m(k) est alors petit par 
rapport a k^, et le vecteur d'onde effectif k cS = k p + i/2£ differe peu du vecteur 
d'onde dans le vide. Par consequent, nous avons 

| 7 | = I - 2lmm{k eS )\ < k 2 ~ k' 2 , \a\, \a\ < I, (3.78) 
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et nous pouvons developper la racine carre dans (13.771) pour trouver : 



rv 2ik' 
Qi ~ 7777 ?2 ~ — (3.79) 
2 k' a 

Le premier terme dans (13.761) qui correspond a la racine qi, represente done une at- 
tenuation exponentielle le long de la direction du vecteur k 

sgngi 6(zlmg 1 ) e [qiZ » sgn 7 6(z 7 ) e _z/£ , - = — = fc , eff pour Q = 0. 

(3.80) 

Le libre parcours moyen qui apparait ici est modifie par rapport au resultat (13.611) parce 
que k' 7^ K, a cause du terme proportionnel a V^m dans (13.701) . 
La deuxieme racine q 2 donne dans cette approximation 

e(az) e mz « e(az) e - |z|/ ^ 2 , £ 2 = 7^7 (3.81) 

elle correspond done a une attenuation tres rapide, sur une echelle en dessous de la 
longueur d'onde dans le milieu. En outre, suivant le signe de la quantite a, l'expres- 
sion (13.811) est non nulle pour des distances z positives ou negatives, done pour des 
sources en amont ou en aval du point d' observation. Dans 1' approximation de BOUR- 
RET, Ton s' attend a 

1 d 2 m 

parce que la partie imaginaire de l'operateur de masse (reliee au libre parcours moyen) 
presente un minimum autour de k fa k (voir la figure l3~Tl) . L expression (13.811) porte 
alors sur les positions z positives (des sources en amont du point d' observation). 

II est justifie de negliger cette deuxieme contribution a la fonction d' attenuation. 
Nous cherchons en effet la portee la plus grande en fonction de la distance du point 
d'observation, et la premiere racine donne une portee beaucoup plus grande que la 
deuxieme. 9 Si nous nous placons a une echelle spatiale plus grande que la longeur 
d'onde, nous pouvons done negliger la contribution de la racine q 2 . Le resultat de 
l'integrale (13.761) sur le vecteur d'onde longitudinal q est alors : 

. expig^ 2iK'sgnz6[z(7 - aQ 2 )} 2 

&q 2k>q- h + iaQ 2 + iaq 2 = 4k" - a( 7 - aQ 2 ) ~ ° Q ^ ] 

(3.82) 

oil nous avons re-exprime 7 en fonction de Q. Nous notons que la fonction assure 
que 1' exponentielle est toujours decroissante en fonction de \z\. 



9 Nous notons aussi que Rytov a calcule la fonction de GREEN moyenne par une technique analogue 
(eq. (4.54) de | Rytov |). II trouve egalement deux contributions, dont Fune s'attenue a l'echelle du libre 
parcours moyen et l'autre a une portee de Fordre de la longueur de correlation £ c . II neglige ensuite 
cette derniere contribution, en se pla5ant a une echelle spatiale plus grande que l c . 
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II reste maintenant a effectuer l'integration sur le vecteur d'onde transverse Q. 
L' integrate prend la forme suivante : 

2ifc'sgnz / W e ^S-"®y*Q* exp[-( 7 - aQ 2 )z/2k'] (3.83) 

II convient de distinguer entre les cas z > et z < (positions sur le rayon lumineux 
anterieures et posterieures au point d'observation). 



Le cas « causal » z > 0. Le module du vecteur d'onde Q est limite a l'intervalle 
< Q < (7/a) 1 / 2 . Nous effectuons l'integration sur Tangle azimuthal du vecteur Q. 
En utilisant une nouvelle variable d'integration, l'integrale s'ecrit 



Mk' 



7/" 

- Jds J (/y/a - sR) e - asz ' 2k \ s = 1 /a- Q 2 (3.84) 


Nous avons egalement neglige le terme proportionnel a aaQ 2 au denominateur 
de (|3.83b . Comme Q 2 est limite par 7/a, ce terme est petit par rapport a k' 2 en vertu 
de l'hypothese (13.781) que nous avons deja utilisee ci-dessus (pour calculer qxp). 

II est facile maintenant de determiner l'echelle caracteristique pour la distance 
transverse R. A partir de 1' argument de la fonction de BESSEL J , nous obtenons : 




6R~J-~ A -2^- ~JU£ (3.85) 



(Nous avons utilise l'echelle de variation \/i c de l'operateur de masse pour estimer 
a ~ £ c (ko/ty/ko = Ic/i.) Avant de donner une interpretation physique de ce resultat, 
etudions le deuxieme cas. 



Le cas « acausal » z < 0. Le nom « acausal » provient du fait, rappelons-le, que les 
positions z < se trouvent en amont du point d'observation : elles decrivent done 
1' influence sur la luminance du milieu que le rayon n'a pas encore traverse d'un point 
de vue geometrique. 

L'integrale (13.831 ) porte maintenant sur les vecteurs d'onde (7/a) 1 / 2 < Q < 00. 
Avec une changement de variables similaire, elle prend la forme 

00 

~8^J dS J o(VV a + SjR ) e~ aslzl/2k ', s = Q 2 - 7 /a (3.86) 


Nous avons encore neglige la contribution en Q 2 dans le denominateur. Sur les dis- 
tances qui contribuent a l'integrale (as\z\/2k r ~ 1), elle est negligeable si z ^> a/2k'. 
Puisque cette distance est bien en-dessous de la longueur d'onde, cette approximation 
est justifiee aux grandes echelles qui nous interessent ici. 
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FIG. 3.2 - Representation de la fonction de phase nonlocale. Les processus de diffusion 
a l'interieur des ellipses a la luminance au point d'observation r. Les ellipses represen- 
ted leur poids respectif par des courbes de niveaux. Elles sont etirees autour du rayon 
arrivant au point r dans la direction d'observation k. La partie « en amont »donne 
egalement une contribution. 

Pour l'echelle transverse 5R, nous trouvons alors deux limites suivant la distance 
\z\ du point d'observation : 



Discussion 

Nous constatons qu'en un point d'observation r donne, la luminance est la somme 
des processus de diffusion qui ont amene de la lumiere dans la direction d'observation 
k et qui ont eu lieu dans un « lobe » autour du rayon geometrique de direction k qui 
aboutit a r. Ce ne sont plus des processus de diffusion qui ont eu lieu exactement sur le 
rayon qui contribuent a la luminance observee : le terme de diffusion de l'equation du 
transfer! radiatif est devenu non-local. La largeur transverse des « lobes »est de l'ordre 
de 5R ~ (4A) 1 / 2 . En outre, une contribution non nulle de la luminance provient de 
positions « acausales » sur le rayon et en aval du point d'observation. La situation est 
illustree sur la figure l3~2l par des courbes de niveaux. 

Nous constatons que la longueur de non-localite transverse est egale a la taille de la 
zone de Fresnel lorsque Ton observe une source ponctuelle de longueur d'onde A a 
la distance £ c de la longueur de correlation. La non-localite de l'ETR que nous mettons 
ici en evidence est done reliee a la diffraction dans le champ intermediate (diffraction 
de Fresnel par rapport a celle de Fraunhofer dans le champ lointain). En outre, 
l'echelle 5R suggere que la deviation par rapport a l'ETR locale est due a la diffusion 
recurrente entre deux diffuseurs qui sont places a l'interieur d'un rayon de correla- 
tion et qui s'eclairent mutuellement par leurs champs proches et intermediaires. C'est 




(3.87) 
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en effet une des predictions de Lagendijk et VAN Tiggelen pour les phenomenes 
physiques au-dela de l'ETR standard [ vanTiggelen96 1 . 

Notre calcul permet d'obtenir encore un autre resultat curieux qui concerne le com- 
portement a grande distance \z\ ^> i de la fonction d' attenuation. Dans cette limite, 
nous pouvons calculer les integrales sur s (I3.84ll3~86l) de facon approchee. Leurs con- 
tributions dominantes proviennent en effet de la region s — > 0, de sorte que nous 
obtenons, pour le cas causal par exemple, 

-i- 7 As U^faR) e— = i^SI (3 .88) 

o 

Si nous pouvons resoudre l'echelle transverse 8R, la fonction d' attenuation est done de 
longue portee, avec une decroissance en 1/z. Nous n'avons pas encore trouve d'inter- 
pretation physique simple de ce resultat. II est d'autant plus surprenant que Ton trouve 
un comportement analogue pour le cas « acausal ». 

II est cependant evident que ce resultat est du aux vecteurs d'onde transverses Q 2 ~ 
7/a, qui s'accompagnent d'une decroissance tres lente de la fonction d'attenuation 
en fonction de \z\ (voir (13.83b "). L'on peut alors revenir a la condition (13.741) sur les 
vecteurs d'onde (Q < \q\ < l/£ c ) et se demander si elle est encore verifiee pour 
la situation que nous envisageons ici. En imposant (7/a) 1 / 2 < l/£ c et en utilisant 
l'estimation pour a donnee ci-dessus, nous trouvons la condition suivante : 

4 3* (3.89) 

Nos conclusions sont done seulement valables lorsque la longueur de correlation est 
plus petite que la longueur d'onde. Ceci n'est pas tres etonnant parce que e'est seule- 
ment dans ce regime que la longueur de non-localite 5R est plus grande que la longueur 
de correlation (sur laquelle il faut toujours moyenner dans une theorie statistique de 
transport). Ce regime apparait frequemment dans le domaine optique lors de la dif- 
fusion par des tres petites particules. II est par contre plus difficile a realiser avec les 
atomes. 

Les atomes se trouvent generalement dans le regime oppose ou la longueur d'onde 
est petite par rapport a la longueur de correlation du milieu aleatoire (regime semi- 
classique). Dans ce regime, l'on n'a pas le droit d'etendre le domaine d'integration a 
des vecteurs d'onde plus grand que \/i c . En limitant les integrales sur Q a cet inter- 
valle, e'est alors la longueur de correlation qui donne l'echelle spatiale minimale pour 
l'ETR : 

regime semi-classique : 5R ~ £ c (3.90) 

En outre, la fonction d'attenuation ne contient qu'une partie causale, avec une portee 
longitudinale donnee par le libre parcours moyen i. 

Conclusion 

L' equation du transfer! radiatif est valable lorsque la luminance varie lentement a 
l'echelle d'une longueur de Fresnel, definie par la longueur d'onde et la longueur de 
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correlation. Les deviations par rapport a l'ETR proviennent alors de la diffusion recur- 
rente de la lumiere par des diffuseurs qui s'eclairement mutuellement par leurs champs 
proche et intermediate. Cette description suppose des diffuseurs petits par rapport a 
la longueur d'onde. Dans la limite opposee (regime semi-classique), la longueur min- 
imale pour l'ETR est donnee par la longueur de correlation. 

Nous avons obtenu ce resultat en calculant une correction a l'equation du transfert 
lorsque Ton se rapproche des echelles spatiales plus courtes. Ce calcul met en outre 
en evidence une correction au libre parcours moyen qui fait intervenir la derivee de 
l'operateur de masse (le passage de k a k' dans (13.701) ). Cette correction est probable- 
ment liee a la vitesse de transport du champ dans le milieu qui n'est donnee ni par 
la vitesse de groupe ni celle de phase [vanTiggelen96|. Elle apparait encore plus net- 
tement dans une approche au transport dependante du temps, que nous n' avons pas 
suivie ici. 

Relation de dispersion non-triviale. Pour finir, une question d'ordre experimental : 
comment est-ce possible que l'operateur de masse presente une variation importante 
avec le vecteur d'onde k ? L'on aurait alors une vitesse de transport k' tres differente 
de celle de propagation libre, ainsi qu'un libre parcours moyen £' fortement corrige. 
On peut esperer que tel est le cas dans les solides ou la relation de dispersion presente 
des fortes deviations par rapport a l'espace libre, notamment aux bords de la zone de 
BRILLUOIN (aux bords des bandes permises). Une autre voie peut etre ouverte par 
des particules diffuseurs avec une « structure interne », des billes dielectriques avec un 
noyau creux, par exemple. En effet, la fonction de correlation de la constante dielec- 
trique g{r) va alors retracer, a courte distance, le profil de l'indice des billes, et ceci se 
traduira par une variation importante de l'operateur de masse avec le vecteur d'onde. 
Ce dernier est en effet donne, dans 1' approximation de Bourret et pour des correla- 
tion isotropes, par 

\k 4 °r 

m (k) = ^ J dr g{r) (e i(fco+fe)r - e^ "^') 
o 

ce qui revient essentiellement a la transformee de Fourier du profil d'indice radial. 
D'un point de vue experimental, l'on pourrait done faire varier la structure interne des 
billes pour observer un comportement different du libre parcours moyen. Malheureuse- 
ment, ceci est une variable qui a egalement une influence sur d'autres quantites (le con- 
traste dielectrique total, moyenne sur la bille, par exemple). Une astuce peut eviter cet 
inconvenient : se placer dans une geometrie en deux dimensions ou le vecteur d'onde 
est fixe par la projection sur le plan de symetrie. En choisisant Tangle d'incidence, l'on 
peut alors faire varier le vecteur d'onde sans changer la frequence. 

II reste a evaluer quel contraste dielectrique et quelle concentration de diffuseurs 
sont necessaires pour que l'effet soit observable. En premiere approche, l'on peut en- 
visager d'utiliser 1' approximation de BOURRET a cet effet. 
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3.5 Reciprocity et retro-diffusion exaltee 



Pour finir, nous voudrions etudier la consequence de la reciprocite de la diffusion 
pour 1' equation de Bethe-Salpeter dans la presente formulation. Sans en donner 
une justification precise, la reciprocite exprime que les deux processus de diffusion 
suivants ont la meme amplitude de diffusion : 

k. — r k et — k — > —n 

parce qu'ils s'obtiennent l'un de l'autre par un simple renversement du temps. Cette 
invariance existe pour notre theorie scalaire, elle n'existe plus pour la diffusion electro- 
magnetique en presence d'un champ magnetique, par exemple. Si les deux processus 
ont la meme amplitude, ceci reste vrai aussi apres la moyenne sur le potentiel aleatoire. 
Par consequent, Ton s' attend a une intensite augmentee (par interference constructive) 
dans ces directions. 

En comparant aux arguments de l'operateur d'intensite dans (13. 501) . nous consta- 
tons qu'il correspond a la correlation entre deux processus reciproques (done de meme 
amplitude) lorsque nous avons 



k+|q = -(rc-iq) 
K + iq = -(k-|q) 



k=-K. (3.91) 



21 V" 2' 

C'est done le cas pour la diffusion vers l'arriere. Dans 1' approximation de BOURRET, 
l'operateur d'intensite prend la forme (13.661) et la probabilite relative de diffusion vers 
l'arriere et vers l'avant est donnee par 

arriere S(2 k) 
avant ~ 5(0) 

(S(k) est la densite spectrale du potentiel aleatoire, supposee isotrope ici.) Cette prob- 
abilite est faible lorsque la longueur d'onde est petite par rapport a la longueur de cor- 
relation. En outre, elle ne presente pas de trace particuliere d'un effet d' interference 
dans la direction arriere. Ceci est du au fait que 1' approximation de BOURRET ne prend 
en compte que les processus de simple diffusion. Par contre, si Ton utilise une expres- 
sion pour l'operateur d'intensite qui aille au-dela de 1' approximation de BOURRET, 
Ton peut tres bien decrire la retro-diffusion exaltee dans le formalisme de l'equation 
du transfer! radiatif. La bonne nouvelle est qu'il suffit d'inclure une serie certes infinie 
de diagrammes (les maximally crossed diagrams), mais que leur calcul necessite pas, 
dans la pratique, d'aller au-dela de 1' approximation de BOURRET : Ton s'en sort done 
en « ameliorant le formalisme ». Nous renvoyons aux travaux de TSANG [Tsang85| 



ainsi qu'au § 6 du livre de Ping Sheng [Sheng95 1 pour plus de details. 



3.6 Conclusion et perspectives 

L'equation du transfert radiatif ne s'en sort pas mal de cette discussion : sa forme 
traduit bien le transport de la luminance (la transformee de Wigner du champ) aux 
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grandes echelles spatiales, meme si le potentiel des diffuseurs est fort et que 1' approx- 
imation de Born echoue. II convient cependant de rappeler que l'ETR generalisee a 
ete obtenue pour un milieu statistiquement homogene, done sans prendre en compte les 
effets de bords ni les interfaces. II faut done modeliser ceux-ci au cas par cas avec des 
coefficients de reflexion et de transmission, par exemple. En outre, l'ETR est seulement 
valable lorsque Ton peut parler de la propagation d'ondes dans le milieu : formelle- 
ment, ceci se traduit par un libre parcours moyen beaucoup plus grand que la longueur 
d'onde (regime de faible diffusion). L'ETR ne donne done aucun acces a la localisation 
des ondes. 

Dans la litterature, Ton rencontre souvent des affirmations moins ambitieuses quant 
au statut de l'ETR : celle-ci ne serait valable que dans 1' approximation de BOUR- 
RET llKravtsov961 ILuck93l ou de "single group". En effet, l'analyse de BARABA- 
NENKOV et Finkel'berg I Barabanenk ov67l montre que la validite de l'approxima- 
tion de BOURRET entraine le regime de faible diffusion : le champ se propage dans 
le milieu, les diffuseurs se voient dans le champ lointain et la dispersion spatiale (la 
dependance de l'operateur de masse avec le vecteur d'onde) est negligeable. Sous ce 
point de vue, 1' approximation de BOURRET apparait done comme une condition suff- 
isante pour l'application de l'ETR. Le libre parcours moyen et la fonction de phase 
sont alors donnes par les operateurs de masse et d'intensite que Ton calcule en ten- 
ant compte d'un groupe de diffuseurs avec k = 1, 2, ... particules (caracterise par une 
fonction de correlation a k points). Cette "single group approximation" revient a 1' ap- 
proximation de Born (les formules de Ryzhik, Papanicolaou et Keller) pour 
la diffusion par un potentiel aleatoire continue. Pour des particules diffuseurs, des cor- 
relations entre particules peuvent done en principe etre decrites par le formalisme, mais 
il semble difficile dans la pratique de calculer les sections efficaces de diffusion (voir 
[ Tsang85{ IMacKintosh89l pour des exemples). 

L'equation du transfert radiatif generalisee (« ETR* », eq. 13.561) que nous avons 
trouvee ici ne depend pas, a priori, de 1' approximation de BOURRET. Elle va done au- 
dela de 1' approximation de la diffusion simple, tout en supposant que la diffusion est 
« faible » (ondes non-localisees). Les travaux des Russes indiquent les effets physiques 
nouveaux qui sont decrits par l'ETR* : 

- la diffusion devient dependante (diffuseurs s'eclairant dans le champ proche et 
intermediaire, diffusion recurrente) 

- la relation de dispersion s'elargit (dispersion spatiale : la partie imaginaire de la 
fonction de Green moyennee n'est plus une fonction 5) 

- il apparait une vitesse de transport pour l'intensite des ondes. 

Nous avons identifie les traces de ces effets en developpant la partie « propagation » de 
l'equation de Bethe-Salpeter au-dela de 1' approximation habituelle des grandes 
echelles. II faut cependant rester prudent sur ces resultats parce que les autres termes 
dans l'equation de B.-S. donnent des contributions analogues, la fonction de Green 
moyennee et l'operateur d'intensite. En outre, ces contributions ne sont pas indepen- 
dantes parce qu'elles sont reliees entre elles par l'identite de Ward (une generalisation 
du theoreme optique), qui est valable a toute echelle. La liste des effets physiques ci- 
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dessus semble quand meme qualitativement correcte ; Lagendijk et van Tiggelen 
dressent en effet un repertoire analogue [ vanTiggelen96 ] . 

Finalement, tout semble indiquer que 1' equation du transfert radiatif ne peut pas 
decrire la localisation (forte) des ondes dans un milieu aleatoire. Notre etude peut 
tout au plus suggerer l'image physique suivante : nous avons constate que l'equation 
de Bethe-Salpeter est non-locale, avec une fonction de phase non-positive qui 
presente des changements de signe. II est alors concevable que dans l'integrale sur 
les processus de diffusion, ceux-ci « interferent de facon destructive » : aux grandes 
distances, le champ ne se propagerait pas et serait localise dans le milieu diffusant. 
Rappelons dans ce contexte que la fonction de Wigner n'est pas forcement positive, 
et que ses « negativites » sont precisement le resultat d' interferences. 
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Conclusion 



Le transport des ondes 

Nous avons etudie dans ce rapport le transport d'une onde (scalaire, electro- 
magnetique ou de matiere) a travers un milieu aleatoire diffusant. L'on peut distinguer 
trois niveaux de description dans la theorie : 

1. le niveau « microscopique », ou Ton tient compte de 1' aspect ondulatoire, en 
etudiant le champ moyen et la fonction de coherence. A ce niveau, les effets 
d' interference sont inclus dans la theorie, dont les equations de base sont celles 
de Dyson et de Bethe-Salpeter. 

2. un niveau « mesoscopique », ou le champ est decrit par une luminance : elle en 
donne l'intensite, distribute en position et selon les directions de propagation. 
La theorie du transfert radiatif, mais aussi la mecanique classique statistique si 
situent a ce niveau theorique. Les equations de base sont celle du transfert radi- 
atif, de BOLTZMANN ou encore l'equation de Fokker-Planck 1 . Le lien a la 
theorie microscopique passe par 1' identification de la luminance avec la fonction 
de correlation (de coherence) du champ en representation de Wigner. L'on ob- 
tient l'equation du transfert radiatif a partir de celle de Bethe-Salpeter dans 
la limite des echelles spatiales grandes par rapport a la longueur de correlation 
du milieu (les diffuseurs se voient en champ lointain). 

3. Finalement, la diffusion dans le milieu redistribue l'intensite de l'onde du 
faisceau collimate incident vers la partie diffuse de la distribution angulaire. 
Aux grandes distances, la distribution angulaire devient isotrope, et le champ 
est decrit seulement par son densite locale d'energie. Sur de telles echelles 
« macroscopiques », le transport de l'onde est gouverne par une equation de la 
diffusion (spatiale). Dans une experience avec les atomes ou le temps d'interac- 
tion est souvent limite, ce sont seulement les atomes les plus lents qui entrent 
dans le regime du transport diffusif. 

Ces differents niveaux theoriques sont represented sur la figure suivante. 

1 II depend du rapport entre la longueur d'onde et la longueur de correlation si c'est l'ETR ou l'equa- 
tion de F.-P qui est bien adaptee au probleme. La deuxieme correspond au regime des petites longueurs 
d'onde ou la diffusion se produit de preference aux petits angles. 
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equation de 
la diffusion 
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luminance isotrope 
A<9 < 1 
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mesoscopique 



diffusion simple 



ETR generalisee 



grande echelle spatiale 



equation de 
Bethe-Salpeter 



microscopique 



Representation schematique des differents niveaux d' approximation pour le transport 

des ondes dans un milieu aleatoire. 



Perspectives 

Les atomes froids places dans les tavelures lumineuses realisent un probleme de 
transport en milieu aleatoire dont on peut facilement varier les parametres caracteris- 
tiques (amplitude du potentiel, longueur de correlation par rapport a la longueur d'onde 
(temperature) des atomes, duree de 1' interaction, poids relatif des effets reactifs et dis- 
sipatifs). Le transport d' atomes represente done une sorte de « laboratoire » ou Ton 
peut etudier des regimes differents du transport. Etant donne la difficulte de refroidir 
des atomes en dessous de la limite de recul (longueur d'onde atomique plus grande 
que la longueur de correlation des tavelures), e'est d'abord le regime mesoscopique du 
transport que Ton explorera. Ce n'est pas un domaine depourvu d'interet : Ton peut 
meme envisager de modeliser, en presence du champ de pesanteur, un modele pour la 
conduction des electrons dans un solide avec des defauts (potentiel aleatoire plus force 
constante). 
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La structure interne des atomes offre une autre possibilite interessante : il est pos- 
sible d'etudier par exemple le transport d'un ensemble de spins polarises a travers 
un champ magnetique aleatoire (statique). Ce probleme est pertinent pour l'inter- 
ferometrie atomique a etats internes : un champ magnetique residuel induit en effet 
une rotation aleatoire du spin et reduit le contraste des franges d' interference. D'un 
point de vue theorique, Ton pourra exploiter l'analogie au transfert radiatif de la lu- 
miere polarisee. Le formalisme de Ryzhik, PAPANICOLAOU et KELLER est egale- 
ment suffisamment general pour trouver rapidement les equations de transport du spin 
atomique dans 1' approximation d'un faible champ residuel. 

Dans le domaine du transfert radiatif proprement dit, les considerations de ce 
rapport indiquent que 1' equation du transfert radiatif est suffisamment generale pour 
decrire le regime de la diffusion dependante, a condition de generaliser la section effi- 
cace de diffusion a la diffusion par des agglomerats de particules. De la meme facon, 
l'ETR peut incorporer dans une certaine mesure des effets d' interference comme la 
retrodiffusion exaltee. Elle se limite cependant a une description a grande echelle spa- 
tiale. A plus courte echelle, le transport est regi par une equation non-locale, et il 
n'est plus possible de definir une relation de dispersion pour les ondes dans le milieu. 
L'ETR est egalement limitee au regime de faible diffusion ou le libre parcours moyen 
est beaucoup plus grand que la longueur d'onde. Le regime oppose correspond a la 
localisation des ondes pour laquelle il faut revenir a une description microscopique a 
l'aide de l'equation de Bethe-Salpeter. II est etonnant qu'il faille en elaborer une 
theorie a un niveau si fondamental bien que la localisation se manifeste par 1' absence 
de transport a tres grande echelle spatiale (voir le titre de l'article d' Anderson). II 
semble que se trouve la une explication pour la difficulte notoire du probleme de la 
localisation qui reste ouvert meme apres quarante ans de recherches ardues. 

Carsten Henkel 
Chatenay-Malabry, juillet 1997 
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